'  ^///^  023 13 


<  j  'j 


AFCRI-63-53 


R£SEARCH  IN  SLECTROMAGNXrnCS. 
FLASMAS  AND  NETWORKS 
by 

Li.B.  Fslten 

Pinal  Repprt  No.  PIBMRI-U12 -b3 
For  Item  I  of 

Contract  No.  Ar-19(604y-4l43 
Project  No.  $.635,  Task  No.  56350 
for 

Electronics  Research  Directorate 
Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
United  States  Air  Force 
L.  G.  Hanscom  Field 
Bedford,  Massachusetts 

February  13j  1963 


i' 


A  S  T  I  .9  . 


i'J. 


!\pf?  )i,\i  pm 


POLYTECHNIC  INSTITUTE  OF  BROOKLYN 
MICROWAVE  RESEARCH  INSTITUTE 

ELECTROPHYSICS  DEPARTMENT 


•  Requests  for  additional  copies  by  Agencies  of  the  Department 
of  Defense,  their  contractors,  and  other  Government  agencies  should 
be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA 
services  or  have  their  'need-torknow'  certified  by  the  cognizant 
military  agency  of  their  project  or  contract. " 

"All  other  persons  and  organizations  should  apply  to  the: 

U.  S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHINGTON  25,  D.  C.  " 


AFCRL-63-53 


Report  No.  PIBMRI-m2.63 
Contract  No.  AF-19(  604) -4143 


FINAL  REPORT 
-  ITEM  I  - 

RESEARCH  IN  ELECTROMAGNETICS, 
PLASMAS  AND  NETWORKS 

by 

L.  B.  Felsen 


Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,  New  York 


Final  Report  No.  PIBMRI-1112-63 
Contract  No.  AF-19(604)-4143 
Project  No.  5635,  Task  No.  56350 

February  13,  1963 


Title  Page 
Acknowledgement 
Table  of  Contents 

54  Pages  of  Text  _ _ 

Submitted  by;  I  >  ^  .. 

■L.B.  Pelsen - 

Professor 


Preparer!  for 

Electronics  Research  Directorate 
Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
United  States  Air  Force 
L.  G.  Hanscom  Field 
Bedford,  Massachusetts 


I 


PlBMJU-1112-63 


Acknowledgement 

The  work  reported  herein  waa  supported  by  the  Electronics  Research 
Directorate  of  the  Air  Force  Cambridge  Research  Laboratories*  Office  of  Aerospace 
Research*  (USAF)  Bedford*  Massaebusetts*  under  Contract  No.  AF-19( 604) -4143. 


I 


PmMRZ-lU2-63 


Table  of  Contents 


Acknowledgement 


A.  Preface  1 

B.  Radiation  and  Diffraction  Problems  in  Anisotropic  Regions  3 

1.  Introduction  3 

2.  Uniaxially  anisotropic  media  7 

a.  General  remarks  7 

b.  Radiation  in  the  presence  of  a  perfect  reflector  8 

c.  Radiation  from  a  half-space  9 

d.  Ring  source  and  plasma  column  10 

e.  Diffraction  by  objects  oriented  perpendicular  to  the  optic  axis  11 

3.  Plasma  under  the  influence  of  a  finite  external  magnetic  field  11 

a.  Network  formulation  11 

b.  The  radiation  fields  in  an  infinite  medium  12 

c.  Radiation  in  the  presence  of  a  plasma  half-space  15 

4.  Modes  in  a  hydromagnetic  waveguide  16 

C.  Propagation  and  Diffraction  Problems  in  Isotropic  Plasmas  19 

1.  Piecewise  homogeneous  regions  19 

a.  The  electromagnetic  field  of  a  source- excited  plasma  layer  19 

1'.  The  electromagnetic  spectrum  of  an  infinite  plasma  slab  19 

2'.  The  near  field  and  radiation  pattern  of  a  slot-excited 

plasma  slab  22 

3'.  Radiation  from  semi-infinite  slot-excited  plasma  sheath 

configurations  25 

b.  Radiation  and  propagation  in  the  presence  of  a  compressible 

plasma  half- space  28 

c.  Surface  waves  on  plasma' transition  layers  29 

2,  Inhomogeneous  media  30 

a.  Diffraction  by  objects  in  a  certain  variable  plasma  medium  30 

b.  Determination  of  plasma  density  from  diffraction  measurements  32 

c.  Reflection  and  transmission  of  electromagnetic  waves  in  the 

presence  of  electron  density  gradients  35 


PIBMR1-U12-63 


Table  of  Contenf  (page  2) 


O.  Guiding  and  Scattering  Propertlee  of  Variable  Impedance  and  36 

Abeorbing  Surfacee 

1.  Introduction  36 

2.  Wedge  or  cone  with  a  linearly  varying  surface  impedance  37 

a.  Formulation  emphasizing  surface  guiding  properties  38 

b.  Formulation  emphasizing  scattering  properties  39 

3.  Circular  cylinder  with  peripherally  varying  surface  impedance  40 

4.  Plane  surface  with  bilinear  impedance  variation  41 

5.  Diffraction  by  an  absorbing  half-plane  42 

E.  Other  Diffraction  Problems  44 

1.  Plane  wave  back- scattering  from  a  semi-infinite  cone  44 

2.  Relation  between  a  class  of  two-dimensional  and  three-dimensional 

diffraction  problems  44 

F.  Cerenkov  Radiation  45 

G.  Network  Theory  47 

H.  Publications  and  Talks  50 


f 

I 


PIBMRI-1112-63 


1 


A.  Preface 

This  Final  Report  summarizes  the  work  performed  during  the  period  from 
July  1,  1958  to  December  31,  1962,  under  the  Electromagnetics  and  Networks  phase 
(Item  I)  of  Contract  No.  AF-19( 604) -4143.  The  other  phases  have  already  been  ter¬ 
minated  and  are  described  in  an  earlier  Final  Report. 

The  scope  of  the  research  activities,  theoretical  in  nature,  can  be  assessed 
from  the  Table  of  Contents.  Summary  descriptions  are  given  in  the  text  and  reference 
is  made  to  the  more  detailed  treatments  which  have  been  reported  elsewhere.  Certain 
research  phases  which  have  not  been  completed  will  be  continued  under  the  successor 
Contract  No.  AF-19(628)-2357.  Two  separate  programs  may  be  dis.tinguished:  electro¬ 
magnetics  and  network  theory.  The  electromagnetics  program  has  been  concerned 
with  a  broad  study  of  wave  propagation  and  diffraction  phenomena  in  layered  or  inhomo¬ 
geneous  rhedia;  in  anisotropic  media;  and  by  variously  shaped  objects  having  constant 
or  variable  surface  impedance  properties.  Many  of  these  analyses  have  been  motivated 
by  current  topics  involving  plasmas,  and  such  practical  applications  as  radio  communi¬ 
cation  with,  or  radar  detection  of,  plasma  sheathed  space  vehicles,  plasma  diagnostics, 
etc. ,  have  been  kept  in  mind  when  devising  idealized  prototype  problems  which  can  be 
subjected  to  rigorous  analysis.  A  variety  of  novel  wave  phenomena  have  been  disdovered 
during  these  investigations  and  are  described  at  the  appropriate  places  in  the  text. 

The  networks  program  summarized  in  this  report  covers  many  different 
areas.  While  some  of  the  results  are  extensions  of  already  well-established  classical 
theorems,  others  contain  answers  to  original  and  previously  unsolved  problems.  In  i 
particular,  a  firm  structure  for  Darlington  synthesis  has  finally  been  erected.  Also, 
the  outstanding  problem  of  optimum  multiple-channel  filter  design  has  been  solved 
completely  for  the  case  of  rational  spectral  densities.  This  latter  result  supersedes 
the  work  in  Wiener's  book  and  in  most  Of  the  periodicals  devoted  to  smoothing  and 
prediction.  Most  topics  are  discussed  briefly  with  emphasis  mainly  on  their  pertinent 
contributions. 

The  principal  contributors  to  Phase  I  of  this  Contract  have  been: 


N.  Marcuvitz,  "Theoretical  and  Exploratory  Experimental  Research  in  Electro- 
maghetics.  Networks  and  Related  Solid  State,  and  Plasma  Topics",  Microw.  Res. 
Inst.,  Pblytech.  Inst,  of  Brooklyn,  Final  Report  PIBMRI-806- 60,  Sept.  I960. 
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B.  Radiation  and  Diffraction  Problems  in  Anisotropic  Regions 
1,  Introduction 

While  problems  of  electromagnetic  plane  wave  propagation  in  anisotropic 
media  have  been  explored  in  some  detail  in  connection  with  light  propagation  through 
crystals,  radio  wave  transmission  in  the  presence  of  the  ionosphere  (gyroelectric 
case)  or  microwave  propagation  through  ferrites  (gyromagnetic  case),  much  less  work 
has  been  done  concerning  the  radiation  from  localized  sources  in  or  near  such  media. 

A  principal  reason  for  this  paucity  of  results  (even  at  the  time  of  initiation  of  this 
Contract)  is  the  complication  of  the  radiation  and  diffraction  processes  encountered 
when  sources  or  obstacles  are  embedded  in  anisotropic  regions  characterized  by  a 
tensor  dielectric  constant  and  (or)  permeability.  Little  incentive  was  present  for 
engaging  in  this  difficult  analysis  since  applications  were  generally  concerned  with 
sources  located  exterior  to,  and  removed  "far"  from,  the  anisotropic  domain.  The 
source  field  could  therefore  be  represented  asymptotically  by  its  local  plane  wave 
constituents,  and  the  latter  were  then  traced  through  the  region  of  anisotropy  (ray¬ 
tracing).  More  recent  applications,  however,  require  the  knowledge  of  the  radiation 
characteristics  from  sources  located  inside  the  medium  or  near  its  boundaries,  and  of 
the  effect  of  perturbing  objects  on  the  radiation  or  diffraction  field.  For  example, 
these  results  are  of  intei^est  for  radio  wave  transmission  to  and  from  satellites  passing 
through  the  ionosphere:  radiation  from  plasma-sheathed  space  vehicles  subjected  to 
an  externally  Impressed  magnetic  field  to  eliminate  communication  blackout;  exci¬ 
tation  of  whistlers;  radiation  from  charges  moving  in  or  near  anisotropic  plasmas  or 
ferrites  (Cerenkov  problem);  plasma  diagnostics.  With  these  and  other  applications 
in  mind,  a  basic  research  program  on  radiation  and  diffraction  in  anisotropic  regions 
was  initiated.  To  render  the  mathematical  analysis  tractable,  the  anisotropic  medium 
was  assumed  to  be  representable  by  a  dielectric  and  (or)  permeability  tensor;  condi¬ 
tions  under  which  actual  plasmas  or  ferrites  can  be  so  represented  are  well-known  in 
1  2 

the  literature  ’  and  need  not  be  stated  here. 

The  simplest  type  of  anisotropy  is  the  uniaxial;  in  an  orthogonal  coordinate 
system,  with  one  of  the  coordinates  parallel  to  the  "optic  axis"  (z-axls)  in  the  medium, 
the  dielectric  tensor  £  and  (or)  the  permeability  tensor  fi  is  diagonal,  with  two  iden¬ 
tical  elements  along  (x,  y)  and  a  different  element  along  z.  In  vector  notation: 


‘lit  +  ‘Z^o^o 


(1) 


where  is  a  unit  dyadic  in  the  domain  transverse  to  z,  and  is  a  unit  vector  in  the 
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positive  z-dlrection.  In  a  plasma  or  ferrite  subjected  to  a  steady  magnetic  field 
H  »  z  H  ,  the  uniaxial  case  obtains  in  the  limit  of  infinite  H  ,  when  the  electrons 
are  constrained  to  move  paralL.i  to  the  z-axis  only.  If  tlie  region  is  plane  stratified 
along  z,  and  transversely  eithc-  n  :houpf<ed,  or  bounded  by  a  perfectly  conducting 
cylindrical  surface  (of  arbitrary  cross  section)  whose  axis  Is  parallel  to  z,  it  is  not 
difficult  to  show  that  the  fields  radiated  by  arbitrary  source  distributions  can  be  de¬ 
composed  into  conventional  E  and  H  modes  relative  to  the  z-axis;  the  modal  amplitudes 
satisfy  transmission  line  equations  which  are  only  slightly  different  from  those  fa¬ 
miliar  from  the  theory  of  isotropic  waveguides.  Radiation  and  diffraction  problems 
in  such  regions  can  therefore  be  analyzed  and  formally  solved  by  slight  modifications 
of  conventional  methods;  the  asymptotic  investigation  of  the  radiation  field,  however, 
reveals  phenomena  which  differ  drastically  from  those  encountered  in  isotropic  media. 
Especially  interesting  is  the  case  t,<  0  or  ^2^  0  which  gives  rise  to  focusing  and 
shadow  effects  that  are  completely  absent  in  isotropic  regions  (Sec.  2b).  If  the  optic 
axis  is  perpendicular  to  the  obstacle  or  boundary  surface,  certain  tvi'c-dimensional 
diffraction  problems  can  still  be  solved  by  relating  them  to  equivalent  isotropic  prob¬ 
lems  (Sec.  2e).  No  such  simple  methods  can  be  employed,  however,  when  the  optic 
axis  is  inclined  arbitrarily  with  respect  to  the  boundary. 

These  studies  of  radiation  and  diffraction  in  uniaxially  anisotropic  regions 
lend  insight  into  some  of  the  effects  to  be  encountered  under  more  general  conditions 
provided  that  the  results  are  phrased  in  a  manner  which  highlights  the  physical  mech¬ 
anism  of  the  process  involved.  At  distances  several  wavelengths  from  the  source, 
the  rigorous  expressions  for  the  fields  can  be  approximated  asymptotically  in  terms 
of  local  plane  wave  contributions  modified  to  account  for  the  anisotropy.  A  most 
significant  modification  is  the  distinction  between  the  direction  of  the  phase  propaga¬ 
tion  vector  (phase  gradient  vector,  £=  K£Q+Cf£^)  and  the  direction  of  the  mean 

energy  flow  vector  (Poynting  or  ray  vector,  S  =  +  z  „S  )  in  a  propagating  plane 

® ^  o  p  o  z 

wave.  (£^  is  the  transverse  vector  coordinate).  Their  relation  can  be  inferred  from 
the  refractive  index  diagram  n(6)  vs.  0  which  gives  the  value  of  the  refractive  index  n 
for  different  observation  angles  6  with  respect  to  the  optic  axis:  the  endpoints  of  the 
vector  £  describe  the  n(0)  vs.  6  surface,  and  the  vector  S  is  perpendicular  to  the  sur¬ 
face.^  In  an  isotropic  medium  (Cj  =  €2.  refractive  index  surface  is  a 

sphere  whence  £  and  S  are  parallel.  In  a  lossless  homogeneous  medium  with  a  uniaxial 
dielectric  tensor  ^  plasma  under  the  influence  of  an  infinite  ex¬ 

ternal  magnetic  field,  is  positive  and  niay  be  positive  or  negative.  In  the  former 
case,  for  E  modes  along  z,  the  n  vs.  0  surface  is  an  ellipsoid  while  in  the  latter  case, 
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{a)*2^0  «2^° 

Fig.  1  -  Refractive  index  curves  for  uniaxial  plasma  medium 


£  and  S  are  no  longer  parallel,  and  when  0,  the  range  of  propagating  rays  is 
restricted  to  the  conical  region  0<  0^  ;  the  region  Tr/2>  e>  0^  is  then  a  shadow  region 
which  has  no  counterpart  in  isotropic  problems.  (These  remarks,  made  for  the  half- 
space  z>0,  apply  also  to  z<0  since  the  refractive  index  surface  is  symmetrical 
about  the  z  -  0  plane). 

The  refractive  index  curves  provide  a  pictorial  representation  of  the  propa¬ 
gation  characteristics  of  the  anisotropic  medium  which  is  expressed  analytically  by 
the  dispersion  relation 

H(k.o)  =  0  ,  (2) 


where  K  and  a  are  the  longitudinal  and  transverse  components,  respectively,  of  the 

2  2  2 

phase  propagation  vector  £.  For  the  uniaxial  medium,  (2)  is  simply  +  (o  /e)  =  , 

e  B  *2^‘l*  obtains  the  simple  surfaces  in  Fig.  1.  For  a  gyrotropic 

plasma  (subjected  to  a  finite  external  magnetic  field  H^),  where  the  dielectric  tensor 


has  the  form 


H  modes  along  z  are  not  affected  by  the  anisotropy,  and  the  H  mode  refractive  index 
surface  is  a  sphere. 
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the  dispersion  relation  (2)  is  much  more  complicated  than  for  the  uniaxial  case,  and 
the  refractive  index  curves  are  separated  into  some  eight  categories^  distinguished 
by  the  magnitude  of  the  ratios  and  •  where  u,  and  are  the  applied, 

plasma,  and  cyclotron  frequencies,  respectively  (in  a  simplest  description  of  a  gyro- 
tropic  plasma,  only  the  electrons  are  considered  mobile  and  collisions  and  compres- 
slonal  effects  are  neglected;  in  this  case,  e.  ,  ,  in  (3)  are  real  functions  of  u,  w  ,  w  ). 

if  5  pc 

However,  the  refractive  index  surfaces  still  characterize  the  relation  between  p  and 
S  as  above  and  therefore  contain  much  pertinent  information  about  the  plane  wave 
propagation  characteristics  in  the  medium.  The  n  vs,  0  plots  also  provide  the  means 
of  determining  the  ray  directions  across  surfaces  of  discontinuity  arising  from  abrupt 
changes  in  the  medium  properties.  ^ 

A  meaningful  interpretation  of  radiation  and  diffraction  phenomena  in  aniso¬ 
tropic  regions  may  therefore  be  achieved  by  expressing  the  radiation  (or  diffraction) 
field  in  a  form  identifiable  in  terms  of  incident,  reflected,  refracted,  diffracted,  etc., 
rays  whose  progress  through  the  medium  is  governed  by  the  plane  wave  refractive 
index  curves.  The  rigorous  solution  of  special  problems  furnishes  the  excitation 
amplitudes  of  the  various  ray  contributions,  and  their  further  progress  can  then  be 
charted  by  the  above-sketched  ray-tracing  procedure.  While  the  ray-tracing  technique 
is  adequate  to  account  for  what  may  be  called  the  geometric-optical  field,  it  does  not 
yield  information  on  interface  phenomena  associated  with  surface  waves,  lateral  waves, 
etc.,  or  on  transition  phenomena  encountered  when  passing  from  one  geometric -optical 
wave  domain  to  the  other.  These  effects  must  likewise  be  assessed  from  a  rigorous 
analysis  of  tractable  prototype  problems.  Thus,  there  is  a  need  for  studying  "canon¬ 
ical"  radiation  and  diffraction  problems  in  anisotropic  media  to  provide  the  type  of 
information  available  for  similar  problems  in  isotropic  regions  -  both  in  the  low- 
frequency  (Rayleigh)  and  high-frequency  (quasi-optic)  domains.  The  difficulties  en¬ 
countered  even  in  isotropic  media  in  the  "resonance  region*  intermediate  between 
the  low-  and  high-frequency  domains  are  expected  to  be  further  compounded  when 
anisotropy  is  present. 

To  gain  a  quantitative  understanding  of  the  problems  outlined  above  as  ap¬ 
plied  to  ionized  plasmas,  a  systematic  research  program  was  initiated  under  this 
contract.  In  line  with  the  philosophy  that  a  thorough  understanding  of  phenomena  in 
relatively  simple  idealized  situations  will  lead  to  a  fruitful  subsequent  study  of  more 
difficult  problems  which  presently  appear  intractable,  the  anisotropic  plasma  em¬ 
ployed  in  the  analysis  has  so  far  been  either  of  the  uniaxial  or  of  the  simple  gyro  - 
tropic  type  as  in  (1)  and  (2),  respectively.  Even  these  problems  involve  considerable 
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nutthematical  eomplexityi  eapecially  when  the  external  magnetic  field  is  inclined 
arbitrarily  with  respect  to  interfaces  or  perturbing  objects,  and  much  work  remains 
to  be  done  In  this  category.  Results  obtained  from  studies  on  this  contract  are  sxun- 
marlzed  below. 

2.  Udiaxiallv  anisotropic  media 
a.  General  remarks 

It  was  pointed  out  in  the  Introduction  that  the  uniaxial  medium  contains  the 
simplest  type  of  anisotropy,  and  that  it  is  possible  to  develop  a  modal  (guided  wave) 
formalism  for  such  regions  by  slight  extensions  of  procedures  employed  in  isotropic 
problems,  provided  that  the  optic  axis  is  parallel  to  the  boundaries  of  any  perfectly 
conducting  cylindrical  object  immersed  in  the  medium.  These  resxilts  were  estab¬ 
lished  in  reference  3,  with  the  additional  generalization  of  longitudinally  variable 
^(z)  and  j|f(z).  It  was  shown  that  problems  involving  non-variable  e  and  ji  are  trans¬ 
formable  directly  into  equivalent  isotropic  problems,  and  that  the  radiation  from 
elementary  source  distributions  (such  as  dipoles,  line  sources,  etc.)  in  an  infinite, 
homogeneous,  uniaxially  anisotropic  space  can  be  expressed  in  closed  form  just  as 
in  an  isotropic  region.  In  the  formulation  of  radiation  problems,  due  regard  must  be 
given  to  a  radiation  condition  which  requires  the  outward  flow  of  energy  from  the 
source  domain,^*' ^  Since  the  directions  of  phase  and  energy  propagation  differ  in 
the  anisotropic  medium,  this  condition  is  not  equivalent  to  the  conventional  one  re¬ 
quiring  an  outwardly  progressing  phase,  and  one  may  encounter  backward  wave  ef¬ 
fects  where  the  pertinent  components  of  the  phase  and  energy  propagation  vectors 
are  oppositely  directed.  It  is  found  that  the  far  field  of  a  lumped  source  comprises 
radially  outgoing  rays  but  that  the  corresponding  phase  fronts  generally  move  in  some 
other  direction  which  can  be  inferred  from  the  refractive  index  plots.  These  state¬ 
ments  remain  true  for  the  gyrotropic  medium. 

An  important  feature  of  the  point  source  radiation  field  in  a  medium  char¬ 
acterized  by  an  open-branched  refractive  index  curve  as  in  Fig.  1(b)  should  be  em¬ 
phasized.  On  the  ray  shadow  boundary  6=0^,  the  fields  are  found  to  diverge;  the 
divergence  is  so  strong  that  the  real  power  flow  through  a  cross-section  surrounding 
the  direction  6  ■  0^  is  also  infinite.  While  this  anomaly  is  removed  if  dissipation  is 
present,  its  resolution  for  the  lossless  case  is  a  matter  of  some  interest  and  was 
actually  studied  first  on  this  contract  in  connection  with  the  more  general  problem  of 
radiation  in  a  plasma  subjected  to  a  finite  external  magnetic  field  (see  (3)  ).^*’  ^  It 
was  found  that  the  above  -  mentioned  singularities  are  due  to  source  distributions 
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whose  strength  drops  to  zero  abruptly,  and  that  the  fields  are  botinded  if  the  source 
intensity  decays  gradually.  Thus,  by  shaping  the  source  distributiotv,  it  is  possible 
to  focus  energy  along  the  direction  of  the  shadow  boundary  cone.  The  radiation  pat¬ 
terns  of  various  source  distributions  are  now  under  study. 

From  the  preceding  remarks,  it  is  noted  that  the  formal  solutions  for 
radiation  from  sources  in  the  presence  of  plane  interfaces  perpendicular  to  z,  or  of 
cylinders,  wedges,  half-planes,  etc. ,  parallel  to  z,  can  be  constructed  with  compara¬ 
tive  ease  for  the  uniaxial  medium.  The  real  task  is  an  asymptotic  evaluation  and 
meaningful  physical  interpretation  of  these  results.  Such  a  study  is  in  progress  and 
has  been  completed  so  far  for  the  case  of  radiation  from  a  plasma  half-space. 

b.  Radiation  in  the  presence  of  a  perfect  reflector 

If  the  refractive  index  surface  has  an  unbounded  branch  as  in  Fig.  1(b),  the 
propagating  rays  due  to  a  longitudinal  electric  current  source  are  confined  to  the 
interior  of  the  cone  6<0^,  where  0  is  measured  from  the  dipole  axis,  and  the  region 
0>  0^  is  opaque  and  supports  only  exponentially  decaying  waves.  The  extent  of  the 
shadow  zone  can  be  reduced  by  insertion  of  a  reflector  (Fig.  2)  which  permits  the 
penetration  of  the  reflected  rays  into  part  of  the  shadow  region  of  the  source.  If  the 

reflector  is  oriented  perpendicular  to  the  z-axis,  its  effect  can  be  rigorously  accounted 
z 


(a)  reflector  perpendicular  to  z-axis  (b)  reflector  inclined  to  z-axis 
Fig.  2  -  Plasma  half-space  bounded  by  perfect  reflector  - 


Rays  from  the  source  penetrate  region  A  directly 
J  but  region  B  only  by  reflection. 

i 

I 
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for  in  terms  of  an  image  source  as  in  Fig.  2(a).  ^  For  an  inclined  reflector,  the  analy¬ 
sis  is  more  complicated  and  has  been  carried  out  so  far  only  for  the  two-dimensional 
problem  of  excitation  by  a  line  source  of  magnetic  currents  directed  transverse  to  z 
(say,  along  the  x-axis).  In  this  case  (Fig.  2(b)),  the  reflected  rays  again  appear  to 
come  from  an  image  source  which  is,  however,  displaced  from  the  mirror  image  point 
by  the  amount  Aj|u'|/A2i  where  A2  =  cos^a  +  (l/£)sin^a,  =  ^1  -  (l/t)j  oin2o, 

^  ~  angle  between  the  reflecting  plane  and  the  y-axis,  and  |u'|  is  the 

distance  from  the  source  to  the  reflector  at  u  =  0.  These  results  are  proved  via  two 
different  methods  in  references  8  and  9. 

c.  Radiation  from  a  half-space 

A  detailed  study  of  radiation  from  a  uniaxially  anisotropic  plasma  half-space 
was  carried  out  for  the  cases  where  the  interface  is  either  perpendicular  to  the  z- 
axis^'  or  inclined  arbitrarily.  For  the  former,  arbitrary  source  distributions  are 
admitted;  the  latter  analysis  has  been  performed  so  far  only  for  the  two-dimensional 
problem  of  magnetic  line  source  excitation.  An  asymptotic  evaluation  of  the  formal 
representation  integrals  by  the  steepest  descent  method  reveals  that  the  far  field  in 
the  interior  or  exterior  of  the  plasma  half-space  comprises  contributions  which  can  be 
identified  as  direct,  reflected  or  refracted  rays  with  the  reflection  and  refraction  laws 
exactly  as  predicted  from  the  refractive  index  curves  (the  reflected  field  can  also  be 
thought  of  as  arising  from  an  image  source  as  in  Fig.  2).  Thus,  it  is  verified  that  these 
fields  can  be  predicted  by  ray-tracing  methods,  once  their  starting  amplitude  is  known. 

However,  there  exists  an  additional  field  contribution  due  to  a  branch  cut 
integral  in  the  asymptotic  formulation  which  can  be  interpreted  as  a  new  type  of  lateral 
wave ;  this  wave  propagates  parallel  to  the  interface  in  the  exterior  region,  sheds  energy 
into  the  plasma  by  critical  refraction  and  is  itself  excited  by  a  ray  incident  from  the 
source  at  the  critical  angle.  For  the  medium  having  a  refractive  index  curve  as  in 
Fig.  1(b),  two  lateral  waves  exist  and  their  ray  paths  are  shown  in  Fig.  3.  The  direc¬ 
tion  of  the  various  rays  can  again  be  predicted  from  Fig.  2(b).  Most  significant  is  the 
fact  that  the  lateral  waves  penetrate  the  shadow  regions  in  the  plasma  which  are  inac¬ 
cessible  to  either  the  direct  or  reflected  rays.  While  the  lateral  wave  fields  in  the 
illuminated  region  are  generally  smaller  than  the  geometric -optical  fields,  they  con¬ 
stitute  the  dominant  contribution  in  the  shadow  zone  in  which  the  geometric -optical 
solutions  are  exponentially  small.  Thus,  lateral  waves  may  be  expected  to  play  a 
more  important  role  on  interfaces  bounding  anisotropic  media  than  on  isotropic  bound¬ 
aries;  their  effect  on  the  latter  has  been  studied  in  some  detail  in  the  literature,  es¬ 
pecially  in  connection  with  pulse  propagation,  ^ 
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Fig.  3  -  Lateral  ray  paths 


The  simple  description  above  fails  in  the  vicinity  of  the  transition  regions 
delimiting  the  various  ray-optical  domains*  and  the  more  complicated  field  calcula- 

g 

tions  required  there  have  likewise  been  carried  out.  In  this  connection,  results  were 
derived  for  the  asymptotic  evaluation  of  integrals  possessing  three  evenly  spaced 
neighboring  saddle  points. 

It  is  apparent  from  the  above  and  should  be  emphasized  that  the  rays,  and 
not  the  phase  gradients,  play  the  significant  role  in  the  physical  interpretation  of  radia¬ 
tion  and  diffraction  phenomena  in  anisotropic  media. 

d.  Ring  source  and  plasma  column 

An  analysis  of  the  radiation  from  a  ring  source  of  magnetic  currents  con¬ 
centric  with  a  solid  or  hollow  plasma  column  inside  a  circular  waveguide  has  been 
carried  out.  The  column  is  subjected  to  either  zero  or  infinite  axial  magnetic  field; 
in  the  latter  Instance,  it  becomes  uniaxially  anisotropic.  The  purpose  of  this  study 
was  to  determine  the  amplitude  of  excitation  of  slow  E-type  surface  waves  which  can 
be  guided  either  by  the  plasma  -  air  interface  or  by  the  plasma  interior;  the  results  are 
of  Interest  for  applications  involving  interaction  with  longitudinally  Injected  electron 
beams.  The  eigenvalue  (transverse  resonance)  problem  was  studied  in  detail  and  nu¬ 
merical  calculations  for  the  possible  surface  waves  of  interest,  as  well  as  their  am¬ 
plitudes  of  excitation,  have  been  presented  in  reference  12. 
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e.  Diffraction  by  objecta  oriented  perpendicular  to  the  optic  axi» 

While  the  remarks  in  Sec.  (a)  have  been  pertinent  for  the  formulation  of 
diffraction  problems  when  a  cylindrical  scatterer  (of  arbitrary  cross  section)  is  ori¬ 
ented  parallel  to  the  optic  (2)  axist  another  interesting  class  involves  perpendicular 
orientation.  For  a  curved  obstacle,  e.  g. .  a  circular  cylinder,  the  direction  of  the 
optic  axis  then  varies  over  the  surface,  thereby  giving  rise  to  effects  not  encountered 
previously.  One  may  show  that  two-dimensional  problems  of  excitation  by  a  line 
source  of  magnetic  currents  parallel  to  the  obstacle  axis  are  transformable  into  equiv¬ 
alent  isotropic  problems,  in  which  one  of  the  coordinates  of  the  given  configuration  is 
'stretched*  according  to  y'  -  ./t' y.  Thus,  a  circular  cylinder  in  the  anisotropic  me¬ 
dium  becomes  equivalent  to  an  elliptic  cylinder  in  an  isotropic  medium;  an  inclined 
plane  surface  transforms  into  a  plane  surface  inclined  at  a  different  angle;  etc.  As  in 
the  previously  mentioned  class  of  problems,  the  major  task  here  is  not  the  construction 
of  formal  solutions  but  rather  their  meaningful  physical  interpretation.  Some  pre¬ 
liminary  results  involving  a  plane  reactive  surface,  and  a  perfectly  conducting  wedge 
have  been  given  in  reference  9. 

3.  Plasma  under  the  influence  of  a  finite  external  magnetic  field 
a.  Network  formulation 

The  simplest  model  of  a  plasma  under  the  influence  of  a  2-directed  finite 
external  magnetic  field  involves  the  dielectric  tensor  (3)  which  is  Hermitean  if  losses 
are  neglected.  The  solution  of  the  source-dependent  Maxwell  field  eqiiations  in  this 
tensor  medium  poses  considerable  difficulties  even  for  the  infinite  homogeneous  space, 
to  say  nothing  of  inhomogeneous  or  bounded  regions,  or  regions  containing  scattering 
objects.'*  A  thorough  study  of  radiation  from  arbitrary  source  distributions  in  trans¬ 
versely  unbounded  anisotropic  media  with  plane  stratification  along  the  .-direction  was 
undertaken.  While  a  general,  abstract  formulation  of  problems  of  this  type  has 

been  available^^,  it  was  found  desirable  to  proceed  in  the  present  case  by  a  generaliza¬ 
tion  of  modal  procedures  familiar  from  the  analysis  of  plane  stratified  isotropic  regions. 
The  Maxwell  field  equations  were  reduced  to  their  transverse  form,  the  transverse  (to 
z)  fields  expressed  in  terms  of  a  continuous  superposition  of  orthogonal  vector  modes 
which  comprise  both  "ordinary"  and  "extraordinary*  constituents  (see  (4)  ),  and  the 
modal  amplitudes  derived  from  the  solutions  of  conventional  transmission  line  equa¬ 
tions  (5).  In  this  manner,  the  field  problem  is  reduced  to  a  modal  network  problem 

. .  ■ '  .  . . 

Simplifications  arise  for  certain  two-dimensional  problems  in  which  the  fields  have 
no  2  variation. 
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involving  source-excited  transmission  lines,  with  the  continuity  requirements  of  the 
transverse  fields  across  a  plane  interface  separating  two  different  media  expressed 
in  terms  of  a  network  which  couples  the  ordinary  and  extraordinary  modes  (see  Fig.  4). 
These  coupling  networks  have  also  been  given  for  plane  boundaries  with  constant  sur¬ 
face  Impedance.  While  the  analysis  is  conceptually  simple,  the  actual  e]q>ressions 
for  the  mode  functions  and  especially  for  the  network  elements  are  quite  involved;  they 
have  been  listed  in  reference  7a. 

The  resulting  representation  for  the  electric  field  has  the  form 

E,(x,y,*)  =  /de  /dt,  rv^(z:€,ti)e^(e.Ti)  +  .  ('*) 

where  and  e  ^  are  the  ordinary  and  extraordinary  mode  eigenvectors,  respectively. 

V  and  V  ,  the  corresponding  mode  amplitudes,  are  determined  from  the  transmission 

O  c 

line  equations 

^  =  1/CZI+  V  ,  -g|-  s  IkYV  +  i  ,  z  *  ,  (5) 

which  hold  separately  for  each  mode  and  in  each  of  the  constituent  media,  v  and  1  are 

voltage  and  current  sources  whose  strengths  are  known  in  terms  of  the  applied  electric 

and  magnetic  current  distribution.  K  =  r[)  is  the  modal  propagation  constant  and 

Z  the  characteristic  impedance,  e  _  _  ,  ,  K_  and  Z_  generally  are  multi- 

y  —  o,  e  o,  e  o,  e  o,  e  ' 

valued  functions  of  a  =  ,  and  their  proper  specification  poses  one  of  the 

problems  in  rendering  the  integral  representation  (4)  unique.  A  similar  formulation 
exists  for  the  transverse  magnetic  field  . 

b.  The  radiation  fields  in  an  infinite  medium 

The  preceding  analysis  permits  the  systematic  construction  of  the  formal 
solutions  for  the  fields  radiated  by  arbitrary,  but  prescribed,  sources  in  a  plasma 
medium  consisting  of  piecewise  constant  layers  along  the  z-dlrectlon.  If  the  region 
extends  to  infinity  along  z,  a  radiation  condition  must  be  imposed  which,  as  noted  on 
p.  5,  requires  an  outward  flow  of  energy.  The  above-mentioned  multivaluedness  of 
the  propagation  constant  K  is  intimately  connected  with  the  radiation  condition,  and 
both  an  analytical  and  a  graphical  procedure  have  been  given  for  its  resolution. 
Because  of  the  exceedingly  complicated  analytical  structure  of  the  integrand  in  (4), 
only  two  problems  have  been  studied  in  detail  at  this  time:  a)  radiation  in  an  infinite, 
homogeneous  medium,  and  b)  radiation  from  a  half-space.  While  in  an  initial  attempt 
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FIG.  4.T>«  PHYSICAL  STRUCTURE  AMO  Tl«  ASSOCIATED  NETWORK 

(LONGITUOINAL  GYROTROPIC  AXIS ) 


PlBMlU-m2-63 


14 


at  reducing  the  Integralsi  a  perturbation  procedure  appropriate  to  weak  external  mag¬ 
netic  field  wae  employed^^  •  the  final  analyeie  wae  carried  out  for  the  general 
caee^*"^  and  the  reeulte  are  aummarised  below. 

The  aeymptotic  evaluation  of  the  fielde  radiated  by  a  dipole  eource  in  the 
infinite  medium  wae  carried  out  by  a  saddle  point  procedure.  Only  real  saddle  points 
which  yield  propagating  wave  contributions  were  taken  into  account;  complex  saddle 
points  give  rise  to  exponentially  decaying  fields.  The  number  of  real  saddle  points  is 
found  to  depend  on  both  the  medium  properties  and  on  the  observation  point  location* 
and  is  predictable  from  the  refractive  index  curves  for  the  medium.  For  a  given  ob¬ 
servation  point  located  at  an  angle  6  in  the  p-z  plane,  the  saddle  points  correspond  to 
those  points  on  the  refractive  index  surface  from  which  the  surface  normals  also  point 
in  the  0  direction.  If  the  n  vs.  6  surface  consists  of  two  separate  branches  as  in  Fig.  5 
for  example,  there  may  be  as  many  as  four  saddle  points,  and  each  contributes  a  propa 
gating  ray  from  the  source  to  the  observation  point.  These  rays  individually  carry 


Fig.  5  -  Graphical  interpretation  of  saddle  point  condition 

energy  in  the  radial  direction,  but  due  to  their  different  phase  velocities,  interference 
results  and  the  total  power  flow  need  not  be  radial.  These  Interference  effects  may 
perhaps  be  helpful  in  connection  with  plasma  diagnostics.  The  ray  amplitudes  decay 
inversely  with  distance  from  the  source  as  in  a  spherical  wave,  and  they  are  also 
found  to  be  proportional  to  ,  where  R^  is  th*  radius  of  curvature  of  the  refractive 
index  curve  at  the  1-th  saddle  point.  If  the  rays  correspond  to  a  point  of  Inflection  T, 
R^  and  two  saddle  points  coalesce.  The  necessary  asymptotic  evaluation  then 
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require*  tlmulteneou*  consideration  of  two  neighboring  saddle  points,  and  the  fields 
are  represented  in  terms  of  Airy  functions  rather  than  exponentials.  also  diverges 
when  the  saddle  point  moves  to  infinity  along  an  uhbounded  branch  of  the  n(8)  curve.  In 
that  case,  the  fields  exhibit  true  singularities  which  can,  however,  be  removed  upon 
replacing  the  delta  function  (point)  source  by  a  source  distribution  of  gradually  di¬ 
minishing  strength  (see  comments  on  p.7  ).  Near  the  singular  directions  where  a  ray 
is  perpendicular  to  an  asymptote  of  the  open-branched  surface,  the  divergent  part  in 
the  integral  representation  for  the  fields  may  be  shown  to  be  the  same  as  for  the  uni¬ 
axial  medium  which  can  be  expressed  in  closed  form;  this  recognition  permits  the  de¬ 
tailed  study  of  the  fields  in  the  vicinity  of  its  singularities. 

The  results  of  this  analysis  have  been  applied  to  the  calculation  of  the 
radiation  patterns  of  a  longitudinal  electric  dipole  source  in  a  gyrotropic  plasma.  All 
possible  ranges  of  and  have  been  considered  and  a  typical  calculation  has 

been  performed  for  each  set  of  parameters  which  gives  rise  to  a  characteristically 
different  refractive  index  diagram.  The  numerical  work  was  carried  out  on  an  elec- 
tronic  computer  and  the  resulting  curves  are  given  in  references  7a,  c.  The  great 
diversity  in  the  radiation  patterns  shows  clearly  the  dependence  of  the  radiation  field 
on  the  medium  parameters  u^,  ,  and  on  the  applied  frequency  w. 

c.  Radiation  in  the  presence  of  a  plasma  half-space 

The  preceding  theory  has  also  been  applied  to  the  analysis  of  the  radiation 
from  a  z-directed  electric  current  element  in  the  interior  or  exterior  of  a  longitudinally 
magnetized  plasma  half-space  bounded  by  the  plane  surface  z  =  0.  The  integrands  in 
(4)  are  now  much  more  complicated  than  for  the  infinite  medium,  and  the  asymptotic 
calculation  of  the  fields  has  been  carried  out  only  for  the  geometric-optical  constituents, 
thereby  ignoring  such  possible  additional  contributions  as  surface  waves  or  lateral 
waves  (arising  from  poles  and  branch  points,  respectively),  which  may  be  important 
near  the  interface  or  in  shadow  regions.  A  complete  ray-optical  interpretation  of  the 

^These  phenomena  do  not  necessarily  occur  in  an  actual  plasma  whose  representation 
by  the  simple  tensor  (3)  is  valid  only  in  the  linear,  small  signal  approximation.  How¬ 
ever,  increased  field  strengths  in  a  plasma  do  most  likely  exist  near  the  singular 
directions  noted  above. 

** These  calculations  were  performed  by  E.  Arbel  in  connection  with  his  doctoral  dis¬ 
sertation.  A  revised  and  condensed  version  of  the  first  part  of  the  thesis  was  intended 
to  be  published  Jointly  by  £.  Arbel  and  L.  B.  Felsen.  Due  to  Arbel' s  departure  to  his 
native  country  of  Israel,  this  revision  was  delayed  and  it  was  decided  to  issue  the 
thesis  in  its  original  form  (reference  7a)  in  order  to  make  its  contents  available.  The 
revised  versions  have  now  been  completed  (references  7b,  c)  and  will  be  issued  shortly. 
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7a. 

asymptotic  solution*  for  the  reflected  and  refracted  fields  has  been  given  and  it 
has  been  shown  how  the  refractive  index  curves  can  be  used  to  establish  the  ray  tra¬ 
jectories.  Under  certain  conditions  one  may  encounter  such  phenomena  as  backward 
refraction  and  conical  refraction. 

The  methods  of  analysis  employed  in  the  anisotropic  medium  studies  on 
this  contract)  and  the  results  achieved,  have  been  compared  in  references  7b,  c  with 
related  work  by  other  authors  which  has  become  available  during  the  period  of  this 
Contract.  Almost  all  of  the  other  investigations  have  proceeded  from  a  three-dimen¬ 
sional  Fourier  integral  approach  rather  than  from  the  two-dimensional  modal  technique 
leading  to  the  double  Fourier  integrals  in  (4).  While  the  triple  Fourier  integral  pro¬ 
cedure  has  a  certain  formal  elegance,  the  derivation  of  asymptotic  results  in  an 
infinite  medium,  or  even  the  formal  analysis  of  plane  stratified  media,  requires  its 
eventual  reduction  to  the  two-dimensional  integral  from  which  is  obtained  at  once  from 
the  modal  (network)  approach.  Especially  in  connection  with  the  plane  stratified 
medium  problem,  the  network  procedure  offers  a  systematization  which  is  lacking  in 
other  methods.  These  advantages  have  recently  been  pointed  out  by  Wu^®  in  his 
analysis  of  radiation  in  the  presence  of  an  anisotropic  plasma  slab. 

4,  Modes  in  a  hydromagnetic  waveguide 

J.  Shmoys 

Modes  in  a  hydromagnetic  waveguide  have  been  investigated  by  Newcomb^^ 
Newcomb  finds  two  sets  of  modes  in  a  waveguide  of  circular  cross-sectioh  containing 
a  low-pressure  perfectly  conducting  fluid  in  an  axial  dc  magnetic  field  at  a  p  'essure 
that  is  much  smaller  than  the  magnetic  field  pressure  h  He  then  finds  how 

these  modes  are  perturbed  (to  first  order)  by  taking  into  account  *Ohm's  law"  in  two 
different  forms  (with  a  gyration  term  and  with  a  conduction  term)  and  by  assuming  the 
pressure  to  be  sufficiently  low  to  make  the  speed  of  sound  small  compared  with  the 
velocity  of  the  magnetohydrodynamic  waves.  The  two  sets  of  modes  are  dalled  trans¬ 
verse  electric  and  principal.  Furthermore,  the  author  remarks:  "There  is  ng 
analogous  to  the  TM  modes  of  an  optical  waveguide,  because  these  would  have  a  longi¬ 
tudinal  electric  field. " 

17 

We  have  shown  that  the  equations  governing  the  behavior  of  the  electro¬ 
magnetic  field  in  an  imperfectly  conducting,  perfectly  tenuous  fluid  can  b4  described 
by  Maxwell's  equations  for  a  uniaxial  medium  whose  dielectric  constant  ih  the  longi¬ 
tudinal  direction  is  (c^  4  o/jw),  where  is  the  dielectric  constant  of  fres  space,  a 
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is  the  conductivity  of  the  fluid,  and  u  the  wave  frequency,  while  the  transverse  dielec¬ 
tric  constant  !•  {  ; — - — jr — )  .  where  B  is  the  applied  dc  magnetic  flux 

V  (j«»/ o)  +  By p  / 

density  and  p  is  the  density  of  the  medium.  The  modes  in  a  waveguide  filled  with  a 
uniaxial  medium  do  fall  into  the  sets  (TE  and  TM),  but  in  this  case,  in  the  limit  of 
infinite  conductivity  the  longitudinal  electric  field  of  the  TKl  modes  tends  to  zero. 

The  transverse  field  dependence  of  the  various  modes  is  identical  with  that  in  a  wave¬ 
guide  filled  with  isotropic  dielectric  (see  also  Sec.  2). 
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C.  Propagation  and  Diffraction  Problenn  in  Isotropic  PUimat 


1,  Piecewlie  homogeneous  reglong 

a.  The  electromagnetic  field  of  >  »ource- excited  plathut  layer 
T.  Tamlr  and  A.  A.  Oliner 

The  nuiin  acope  of  thla  research  effort  was  the  investigation  of  the  electro¬ 
magnetic  properties  of  a  planar  homogeneous,  isotropic  plasma  slab,  including  prob  - 
lems  of  radiation  through  plasma  layers  from  finite  sources.  The  purpose  of  this 
study  was  to  clarify  via  this  idealized  model  the  importance  and  influence  of  plasma 
sheaths  which  surround  a  missile  upon  its  re-entry  in  the  upper  atmosphere,  and  thus 
to  obtain  a  proper  understanding  of  the  effect  that  such  an  ionized  layer  will  have  on 
the  communication  to  and  from  the  re-entering  vehicle. 

The  major  aspects  of  this  investigation  may  be  classified  into  the  following 
three  topics; 

1.  The  electromagnetic  spectrum  of  an  infinite  plasma  slab. 

2.  The  near  field  and  radiation  pattern  of  a  slot-excited  plasma  slab. 

3.  Radiation  from  semi-infinite  slot-excited  plasma  sheath  configurations. 
The  principal  features  of  each  of  these  studies  are  described  separately  below. 

1*.  The  electromagnetic  spectrum  of  an  infinite  plasma  slab 
12  4 

This  part  of  the  study  ’  ’  consists  of  an  extensive  investigation  of  all 
the  types  of  waves  which  may  be  supported  by  a  lossless  and  homogeneous  plasma 
layer.  The  geometry  of  the  configuration  is  shown  in  Fig.  6;  the  plasma  slab  is  charac¬ 
terized  by  a  relative  dielectric  constant: 

m 

where  is  the  plasma  frequency  and  w  is  the  frequency  of  the  applied  field. 

z 


Fig.  6  •  Geometry  of  the  homogeneous  isotropic  plasma  slab 
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Ap 


Fig*  7  -  Diepersion  curves  for  surface  waves  along  the  plasma  slab; 

\/ \  Is  shown  versus  either  t  or  X/ \  for  both  the  short- 
8  P  P 

circuit  (sb)  and  open-circuit  (ob)  bisection  cases,  for  differ¬ 
ent  values  of  a  *  Z-ird/ X  ,  with: 

P 

X:  wavelength  in  free  space, 

X^;  guided  wavelength  along  the  plasma  slab, 

X  «  2ir  c/w  , 

P  P 

c:  velocity  of  light  in  free  space. 

The  dotted  line  shows  the  boundary  between  the  forward  wave 
and  the  backward  wave  regions. 
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In  order  to  find  any  poealble  diacrete  contributions  to  the  spectrunii  the 
problem  was  analysed  by  assuming  either  an  electric  or  a  magnetic  plate  at  the  plane 
of  symmetry  z  a  0;  these  solutions  are  then  referred  to,  respectively,  as  short* 
circuit  (sb)  or  open«circuit  (ob)  bisection  solutions. 

It  was  then  found  that  E  mode  surface  waves  may  be  present  along  the 
slab  whenever  is  negative,  as  shown  by  the  dispersion  curves  of  Pig.  7.  The 
field  of  these  surface  waves  decays  away  from  the  plasma-air  interface  in  both  the 
plasma  and  air  regions.  An  interesting  feature  of  these  waves  is  the  fact  that,  for 
some  values  of  the  parameters  involved,  the  waves  may  be  of  the  backward  type,  as 
indicated  in  Fig.  7.  This  feature  is  novel  for  isotropic  configurations  and  it  is  noted 
that  the  surface  waves  discussed  here  are  independent  modes  and  not  backward  space- 
harmonics  such  as  occur  in  periodic  structures. 

In  addition  to  these  surface  waves,  other  modes  with  complex  (rather  than 
pure  real)  wavenumbers  were  shown  to  be  present  for  E  modes  and  for  negative  values 
oftp.  These  waves,  unlike  improper  (non-modal)  leaky  waves,  possess  fields  which 
are  defined  everywhere  in  space  and  decay  in  all  directions  away  from  the  source. 

The  field  contours  for  this  type  of  wave  are  shown  in  Fig.  8,  which  refers  to  the  case 
of  a  short-circuit  (sb)  bisected  slab.  Although,  as  indicated  in  this  figure,  power  flow 
seems  to  be  associated  with  such  a  complex  wave,  actually  two  such  waves  must  be 
present  simultaneously  in  lossless  media  and  they  have  their  power  flow  in  opposite 
directions;  the  net  result  of  the  existence  of  both  waves  is  then  a  cancellation  of  the 
total  real  power  flow,  with  only  reactive  power  being  left  in  the  form  of  stored  energy 
in  the  neighborhood  of  the  source. 

z 


Fig.  8  -  Field  contours  for  a  complex  spectral  wave 
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When  E  modes  with  or  M  modes  (with  any  value  of  t^)  were  considered, 
it  was  shown  that  no  surface  or  spectral  complen  waves  can  exist.  Instead,  by  using 
a  non* spectral  representation,  it  was  shown  that  improper  leaky  waves  may  be  excited; 
these  leaky  waves  have  a  strong  influence  on  the  shape  of  the  radiation  pattern  due  to  a 
source  located  in  or  near  the  slab.  This  aspect  is  discussed  in  the  next  section. 

2*.  The  near  field  and  radiation  pattern  of  a  slot»excited  plasma  slab 

The  information  obtained  from  the  study  described  above  was  used^'  ^ 
to  find  the  radiation  pattern  of  a  plasma  layer  backed  by  a  metal  plate  when  the  excita¬ 
tion  is  in  the  form  of  a  slot  located  in  the  plate,  as  shown  in  Fig.  9,  where  the  slot  is 
represented  by  a  magnetic  line  source.  This  particular  geometry  was  chosen  since  it 
corresponds  to  a  first  Idealized  model  of  an  antenna  located  on  the  wall  of  a  missile 
when  the  latter  is  enclosed  by  a  plasma  sheath. 


t 


Fig.  9  -  Geometry  of  plasma  layer 

It  was  then  shown  that,  when  t^>  0,  the  near  field  is  dominated  by  the 

P 

presence  of  a  single  leaky  wave.  The  air -plasma  interface  was  then  chosen  as  an 
aperture  surface,  atid  the  leaky  wave  along  this  surface  taken  as  the  aperture  distribu¬ 
tion,  in  order  to  compute  the  associated  radiation  pattern.  The  resulting  pattern 
exhibited  a  sharp  major  peak,  as  expected  from  previous  calctilations  of  such  patterns 
obtained  by  steepest  descent  techniques  (in  the  absence  of  the  plasma  layer,  the  radia¬ 
tion  pattern  is,  of  course,  isotropic  with  respect  to  the  polar  angle  0).  This  major 
peak  is  present  at  an  angle  6^  given  closely  by: 


(2) 
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When  X  model  with  or  H  model  (with  any  value  of  t^)  were  coniidered, 
it  wai  ihown  thnt  no  lurfaee  or  ipeetral  complex  wavei  can  exiit.  Znitead,  by  uiing 
a  nen«ipectral  repreientationi  it  wai  ihown  that  improper  leaky  wavei  may  be  excited; 
theie  leaky  wavei  have  a  itrong  influence  on  the  ihape  of  the  radiation  pattern  due  to  a 
louroe  located  in  or  near  the  ilab.  Thii  aipect  ii  diieuiied  in  the  next  lection. 

2',  The  near  field  and  radiation  pattern  of  a  ilot^exclted  plaima  ilab 

3  6 

The  information  obtained  from  the  etudy  described  above  wai  uied  ' 
to  find  the  radiation  pattern  of  a  plasma  layer  backed  by  a  metal  plate  when  the  excita¬ 
tion  ii  in  the  form  of  a  slot  located  in  the  plate,  ai  shown  in  Pig.  9»  where  the  slot  is 
represented  by  a  magnetic  line  source.  This  particular  geometry  was  chosen  since  it 
eorreipondi  to  a  first  Idealized  model  of  an  antenna  located  on  the  wall  of  a  missile 
when  the  latter  is  enclosed  by  a  plasma  sheath. 


t 


Fig.  9  -  Geometry  of  plasma  layer 

It  was  then  shown  that,  when  c^>  0,  the  near  field  is  dominated  by  the 
)>reieaee  of  a  single  leaky  wave.  The  air -plasma  Interface  was  then  chosen  as  an 
aperture  surface,  and  the  leaky  wave  along  this  surface  taken  as  the  aperture  distribu¬ 
tion,  in  order  to  compute  the  associated  radiation  pattern.  The  resulting  pattern 
exhibited  a  sharp  major  peak,  as  expected  from  previous  calculations  of  such  patterns 
obtained  by  steepest  descent  techniques  (in  the  absence  of  the  plasma  layer,  the  radia¬ 
tion  pattern  is,  of  course,  isotropic  with  respect  to  the  polar  angle  9).  This  major 
peak  is  present  at  an  angle  6^  given  closely  by: 

sine^  -  • 


(2) 


PZBMIU-U12-63 


23 


Thit  angle  6^  ie  the  reeult  which  one  would  ei^ect  from  geometric -optical  argumenta; 
however,  it  also  correeponde  to  the  radiation  from  the  leaky  wave  which  ia  moat 
atrongly  excited.  In  fact,  it  may  be  ahown  that,  in  moat  caaea,  the  radiation  pattern 
can  be  found  by  aaaximlng  the  field  at  the  plaama-air  Interface  to  conaiat  only  of  the 
contribution  due  to  thla  leaky  wave.  Thia  ia  exemplified  by  the  radiation  patterna 
ahown  in  Fig.  10  where  la  the  aaymptotically  rigoroua  (ateepeat  deacent)  reault 
for  the  radiation  pattern,  while  Rp(d)  ia  the  reault  obtained  when  one  aaaumea  the 
near  field  due  only  to  the  major  contributing  leaky  wave.  One  then  obaervea  that  other 
field  contrlbutlona  (from  the  apace  wave  and/or  minor  leaky  wavea)  yield  only  email 
correctiona  to  R  (d). 


Fig.  10  -  Radiation  pattern  for  a  ■  1. 0  and  -  0.  l(d/X  ■  0. 168).  The  firat  leaky 

wave  la  atrongly  dominant;  the  other  leaky  wavea  yield  broadaide  radiation 
only. 
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I 


When  the  alnb  thieknaaa  d  it  relatively  large,  other  minor  leaky  waves 
yield  additional  smaller  peaks,  as  shown  in  Fig.  11.  In  that  easst  the  locations  of  the 
peaks  may  also  he  predicted  and  are  associated  with  the  presence  of  these  latter  peaks, 
as  indicated  in  the  figure.  It  then  follows  that,  in  general,  the  features  of  the  radia* 
don  pattern  can  be  directly  related  to  and  easily  predicted  from  the  values  of  the  wave- 
numbers  of  the  leaky  waves.  When  one  calculates  these  wavenumbers  as  funedons  of 
the  frequency  of  operadon,  the  thickness  of  the  slab  or  the  electron  concentradon  in 
the  plasma,  one  can  easily  find  how  some  of  the  major  features  of  the  radiadon  pattern 
change  when  any  one  of  these  parameters  varies. 


Vlg.  11  -  Radiadon  pattern  for  a  «  10. 0  and  >  0. 5  (d^  ■  !•  The  arrows  show 
the  angles  of  maximum  radiadon  for  the  first  four  leaky  waves;  all  the  other 
leaky  waves  contribute  at  broadside  only. 

When  Cp<  0i  the  Held  la  the  plasma  medium  decays  exponentially  so  that  very  little 
radladoB  is  present.  Then,  energy  may  be  carried  along  the  slab  due  to  a  surface 
wave  and  a  small  amount  of  broadside  radiation  is  present;  the  latter  may  be  shown  to 
be  associated  with  the  presence  of  the  spectral  complex  poles  mendoned  in  See.  A. 
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3'.  Radiation  from  ■eml-itiflnite  ■lot»«tcited  pUstha  iheath  coaflguratloin 

81ae«  »  ailtail*  i«  In  fact  finite  and  not  infinite  in  extenti  all  improvement 
over  the  model  eoaeidered  in  8ec.  B  ie  obtained  by  taking  into  account  the  influence  ot 
fLnltenece  or  of  diaeontinuitiea  which  may  be  present.^  In  particular,  the  two  con¬ 
figurations  shown  la  ITlg.  12  were  considered. 

In  the  first  case,  fig.  12(a),  the  plasma  slhb  is  terminated  at  a  plane  AA*  to 
approximate  the  efiect  of  the  electron  coneentratloa  thinning  out  to  negligible  amounts, 
la  the  second  ease,  Fig.  12(b),  the  presence  of  the  end  wall  of  the  missile  is  accounted 
for  by  letting  the  plasma  layer  eontintte  and  fill  the  quarter-space  vacated  by  the  metal. 

8ince  the  near  field  is  dominated  by  a  single  leaky  wave.  It  was  assumed 
that  the  field  at  the  discontinuity  plane  AA*  is  given  by  this  leaky  wave  alone,  for  both 
cases  shown  in  Fig.  12.  The  radiation  field  was  then  calculated  by  means  of  a  Klrchhoff 
integration  over  this  field,  and  the  higher  mode  efiects  at  the  discontinuity  were  dis¬ 
regarded.  The  result  thus  obtained  was  evaluated  numerically  for  many  values  of  the 
parameters  involved.  In  all  these  eases,  it  was  found  that  the  major  radiation  peak 
discussed  in  See.  B  1»  present  for  the  semi-infinite  case  as  well;  other  peaks  may 
then  also  appear  but  they  are  small  in  amplitude  in  most  eases.  In  addition,  some 
end-fire  (8  «  ir/2)  radiation  is  now  present. 

hiost  notably,  however,  it  was  found  that  the  dlfierence  between  the  semi¬ 
infinite  and  infinite  cases  was  strongest  when  the  source  was  located  close  to  the  dis¬ 
continuity  plane  AA'.  When  the  separation  between  the  source  and  plane  AA'  is  only 
one  or  two  wavelengths,  the  pattern  is  strongly  modified!  when  this  separation  is 
greater  than  10  or  20  wavelengths,  the  pattern  is  indistinguishable  from  that  for  the 
Infinite  case. 

These  results  therefore  show  that,  except  for  some  minor  changes  and 
eneept  when  the  source  is  very  near  to  the  discontinuity,  the  radlafion  patterns  for 
these  cases  are  not  appreciably  different  from  those  for  the  Infinite  case.  Ode  there¬ 
fore  concludes  fiuit  the  major  peak  due  to  the  most  strongly  excited  leaky  wave  is  a 
feature  which  dominates  the  radiation  pattern  and,  as  such,  is  not  much  affected  by 
the  presence  of  discontinuities  as  long  as  these  are  not  located  too  close  to  the  source. 


PtBMIU-m2-6l 


27 


Reference* 


1.  A.  A.  Oliner  end  T.  Temlri  *Beckwerd  Wave*  on  leotropic  Pleeme  Slebe", 
Sleetr<mh)ra.  Group  Kiemo.  66,  PlBldRI>930«6lt  Mierow.  Rea.  Inst.,  Polytech. 
Inst,  of  Brooklyn,  1961. 

2.  T,  Tendr,  *The  Blcetromegnetic  Field  of  a  Source*£xclted  Plasma  Slab*, 
!IPIBMRN932«61,  Ktlcrow.  Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Sept.  1961. 

3.  T.  Tamir  and  A.  A.  Oliner,  *The  Influence  of  Complex  Waves  on  the  Radiation 
Field  of  a  81ot-£xcited  Plasma  Slab*,  PIBMRI-959*61,  Mierow.  Res.  Inst., 
Polytech.  Inst,  of  Brooklyn,  Oct.  1961. 

4.  T.  Tamir  and  A.  A.  Oliner,  *The  Spectrum  of  Electromagnetic  Waves  Guided  by 
a  Plasma  Layer*,  P1BMRI*970>61,  Mierow.  Res.  Inst. ,  Polytech.  Inst,  of 
Brooklyn,  Dec.  1961. 

5.  A.  A.  Oliner  and  T.  Tamir,  *Radiation  from  Semi*Infinite  Slot>£xcited  Plasma 
Sheath  Configurations*,  PIBMRI‘>1029  *  62,  Mierow,  Res.  Inst.,  Polytech.  Inst, 
of  Brooklyn,  May  1962. 

6.  T.  Tamir  and  A.  A.  Oliner,  *The  Near  Field  Supported  by  a  Source>Excited 
Plasma  Slab*,  tJRSI>IRE,  May  1961,  Washington,  D.  C. 


PIBMRI-m2-63 


28 


b.  Radiation  and  propagation  in  the  presence  of  a  compressible  plasina 

half-apace 
J.  Shmoys 

Most  of  the  existing  work  dealing  with  the  propagation  of  electromagnetic 
waves  on  plasma  structures  is  based  on  the  assumption  that  only  electrons  are  mobile 
and  that  the  electron  temperature  is  zero.  This  investigation  deals  explicitly  with  the 
effect  of  finite  electron  temperature.  At  a  finite  temperature  an  electron  gas  can 
propagate  longitudinal  waves  with  a  finite  velocity.  If  it  were  possible  to  affect  the 
radiation  of  electromagnetic  waves  to  an  appreciable  extent  by  changing  plasma  tern* 
perature,  then  the  manner  in  which  this  influence  occurs  would  be  of  interest  both 
from  the  point  of  view  of  the  performance  of  electromagnetic  devices  in  the  presence 
of  warm  plasmas  and  from  the  point  of  view  of  possible  diagnostic  applications. 


The  structure  investigated  here  was  a  semi-infinite  plasma  bounded  from 
the  rest  of  space  by  a  stationary  plane  interface.  A  magnetic  line  current  parallel  to 
the  interface,  located  in  free  space  at  a  finite  distance  from  the  boundary,  was  found 
to  excite  plasma  waves  as  well  as  a  surface  wave.  This  surface  wave  also  exists  in 
the  same  configuration  having  zero  electron  temperature,  and  the  electron  temperature 
affects  its  properties  to  some  extent.  Early  in  this  investigation^  it  was  erroneously 
concluded  that  as  a  result  of  the  finiteness  of  temperature,  a  second  surface  wave 
could  propagate;  it  was  later  found  that  this  second  root  of  the  resonance  equation  does 
not  lie  on  the  proper  sheet  of  the  appropriate  Riemann  surface,  and  thus  does  not  cor¬ 
respond  to  a  surface  wave.  At  zero  temperature  the  surface  wave  can  propagate  only 
at  frequencies  which  are  smaller  than  the  plasma  frequency  by  a  factor  of  at  least 
i,  e, ,  w<Wp/  •/z”.  The  effect  of  finite  temperature  is  to  extend  the  pass -band  of  the 
surface  wave  to  infinite  frequencies.  At  low  frequencies,  below  the  zero  temperature 
cut-off  frequency,  temperature  has  little  effect  on  the  propagation  of  the  surface  wave. 
At  higher  frequencies,  the  dispersion  curve  becomes  very  steep  so  that  the  surface 
wave  will  hardly  be  excited.  Almost  all  the  power  carried  by  the  surface  wave  is  in 
the  electromagnetic  field.  A  calculation  of  the  total  power  transfer  into  the  longitudinal 
oscillations  is  still  in  progress. 


2  3 

The  results  of  this  investigation  were  issued  in  a  report  and  published  . 

1.  9th  and  12th  Quarterly  Reports. 

2.  A.  Hessel,  N,  Marcuvitz  and  J.  Shmoys,  "Scattering  and  Guided  Waves  at  an  Inter¬ 
face  Between  Air  and  a  Compressible  Plasma",  PIBMRI-953-61,  Microw.  Res. 

Inst. ,  Polytech.  Inst,  of  Brooklyn,  Sept.  1961. 

3.  A.  Hessel,  N.  Marcuvitz  and  J.  Shmoys,  "Scattering  and  Guided  Waves  at  an  Inter¬ 
face  Between  Air  and  a  Compressible  Plasma",  IRE  Trans.  PGAP,  Vol,  AP-10, 
pp.  48-54,  Jan.  1962. 
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c.  Surface  w»vea  on  traiwltion  layeri 

J.  Shmoyt 

An  interface  between  air  and  plaema  (regarded  aa  a  dielectric  with  negative 
dielectric  conetant)  is  capable  of  aupportlug  (E  mode)  surface  waves.  Since  it  is  im¬ 
possible  to  achieve  a  sharp  boundary  between  the  two  regions,  it  is  of  interest  to  in¬ 
vestigate  the  problem  of  propagation  of  surface  waves  when  the  transition  from  air  to 
plasma  is  gradual.  The  E  mode  problem  for  this  case  of  smooth  variation  is  rather 
difficult;  therefore,  it  was  felt  that  some  Insight  might  be  obtained  from  considering 
semi-infinite  regions  of  air  and  plasma  separated  by  a  homogeneous  slab  of  interme¬ 
diate  density.  This  problem  was  solved  by  the  application  of  the  transverse  reso¬ 
nance  procedure.  It  turns  out  that  the  dispersion  curve  for  the  case  of  an  intermediate 


FREE  SPACE 


SURFACE  WAVE 


Fig.  13  -  Transition  layer  geometry 

density  slab  differs  considerably  from  that  for  the  abrupt  transition.  In  the  latter  case, 
there  is  a  single  low  frequency  pass -band  for  the  surface  wave,  extending  up  to 
w  A/s'  (w- the  plasma  frequency  of  the  semi-infinite  plasma  region).  The  transl- 
tim  layer  has  the  effect  of  shifting  the  cut-off  frequency  of  that  band  downward  and  of 
introducing  another  pass -band  at  somewhat  higher  frequencies  (but  still  below  u  ).  In 
additloa,  if  the  transition  layer  is  thin,  two  distinct  modes  of  propagation  may  exist 
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tM»V  thi  ^I^equeiiey  ot  tlie  bafadt  one  ot  theta  belug  a  backwai^d  wave.  The 

detailed  aualyala  oi  the  tti^ohlem  and  numefleal  Veaulte  wei‘e  iaeued  la  an  £lectV6phyalea 
hlelnoifandttmi  1*18X011-1093^63.^  Typical  diapeiraloa  cuifvei  are  eketched  on  l*lg.  14. 
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IHg.  14  «  Dlepereion  curvee  tot  a  transition  layer 
l>  Inheiiiogeneoai  media 

ai  Oiffractlon  by  objects  In  a  certain  variable  plasma  medlnm 

A  transformation,  also  mentioned  in  Section  82,  can  be  employed  to  relate 
certain  three-dimensional  problems  of  diffraction  by  rotatiohally  symmetric  objects 
In  a  homogeneom  medium  to  equivalent  two-dimensional  diffraction  problems  in  an 
Inhomogeneotte  medium.  The  two  related  configurations  are  schematised  in  Fig.  15. 
The  three-dimensional  problem  in  a  medium  characterised  by  the  constant  wave- 


1.  Pt  kiriek  end  J.  siun^s,  'Surface  Waves  on  Plasma  Slabs*,  iClectrophysics 
MemoraAduim  I5|  Ptfllllll-1093-62,  hlicrow.  Res.  Inst. ,  Folytech.  Inst,  of  Brooklyn, 
Oetober  111  19^2. 
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tmmbcr  involvaa  excitation  of  the  rotationally  eymmetric  acatterer  by  a  ring 
•ouree,  while  the  two>dlmenaional  problem  with  medium  parameter  k(y)  haa  a  line 
aouree  parallel  to  the  obataele  axia.  The  intereating  feature  of  the  latter  claaa  of 
diffraction  problema  ia  firatly,  that  the  frequency  characteriatica  of  the  medium  are 
thoae  of  a  loaaleaa  plaama  with  a  variable  electron  denaity  and  aecondly,  that  the  (y) 
direction  along  which  the  medium  variea  need  not  bear  any  relation  to  the  aurface  of 


•INFINITE 

PLANE 


(a)  THREE'OMENSIONAL  CONFIOUPATION  (b)  TWO-DIMENSIONAL  CONFIGURATION 


Fig.  15  >  Related  diffraction  problema 

the  acatterer.  Thuat  a  atudy  of  the  rigoroua  aolutiona  for  varioui  obataele  configura- 
tiona  permita  the  inveatigation  of  wave  perturbationa  cauaed  by  refractive  index  varia- 
tiona  on  the  acatterer  aurface,  thereby  providing  a  more  general  claaa  of  canonical 
problema  than  thoae  available  in  the  literature. 

The  analyaia  and  aome  preliminary  reaulta  have  been  given  in  reference  1. 


1.  L.  B.  Felaen,  'Oiffractioa  by  Objecta  in  a  Certain  Variable  Plaama  Medium*, 

Xleetrophya.  Oroup  Memo.  S3,  R«7S5>59,  PlB-693,  Microw.  Rea.  Inat. ,  Polytech. 
Inet.  of  Brooklyn,  July  1939. 
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b.  Determination  of  plaama  denaity  from  diffraction  measurementa* 

J.  Shmoys 

The  etandard  method  of  microwave  probing  of  plaemas  is  usually  a  slight 
modification  of  the  diagnostic  method  developed  for  ionospheric  sounding.  While  in 
ionospheric  sounding  the  time  delay  in  reflection  is  measured  as  a  function  of  frequency* 
in  plasma  diagnostics  the  phase  shift  in  either  reflection  or  transmission  is  determined. 
Under  some  restrictions  one  can  reconstruct  the  profile  of  the  medium  from  the  varia¬ 
tion  of  these  quantities  with  frequency.  A  plane  stratified  medium  and  an  incident 
plane  wave  is  usually  assumed.  The  two  difficulties  mentioned  above  are;  (1)  fre¬ 
quency  nriust  be  varied  and  (2)  the  plasma  must  be  plane  stratified.  -The  variation  of 
frequency  over  a  wide  range  is  an  experimental  nuisance  in  the  microwave  range.  A 
planar  configuration  in  a  plasma  is  sometimes  difficult  to  obtain.  Many  plasmas 

which  have  cylindrical  symmetry,  would  have  to  be  extremely  large  so  that  cylindrical 
stratification  could  be  considered  approximately  plane. 

It  is  in  order  to  avoid  these  difficulties  that  the  present  diagnostic  method 
is  being  proposed.  We  assume  a  cylindrical  geometry  and  an  incident  plane  wave  of 
fixed  frequency.  If  we  now  measure  the  differential  scattering  cross-section,  i.  e. , 
the  scattered  wave  amplitude  vs.  angle,  we  can  reconstruct  the  density  profile.  The 
assumptions  are  essentially  the  same  as  in  ionospheric  sounding;  the  variation  of 
electron  density  with  radius  must  be  slow  (fractional  charge  per  wavelength  must  be 
small)  so  that  the  ray  theory  of  propagation  of  electromagnetic  waves  be  applicable; 
the  density  variation  must  be  monotonic;  the  collision  frequency  must  be  sufficiently 
small  so  that  its  effects  are  negligible.  Under  these  assumptions  the  calculation  of 
density  from  the  scattering  data  can  be  broken  up  into  two  parts:  the  calculation  of 
the  ray  geometry  from  the  scattering  data,  and  the  calculation  of  the  electron  density 
distribution  from  the  ray  geometry.  The  relation  between  the  deflection  angle  of  a 
ray  as  a  fxmctlon  of  distance  between  that  ray  and  the  central  ray  (impact  parameter) 
in  the  incident  wave  (cf.  Fig.  16)  is  obtained  from  the  differential  scattering  cross- 
section  by  a  consideration  of  conservation  of  energy. 

We  have  to  distinguish  between  the  somewhat  simpler  case  of  an  impene¬ 
trable  (overdense)  plasma  in  which  case  the  ray  deflection  angle  is  a  monotonic  func¬ 
tion  of  the  Impact  parameter,  and  the  more  involved  one  of  a  penetrable  plasma, 

"'This  topic  was  mentioned  in  the  previous  Final  Report  dealing  with  Contract  Items 
2-4.  It  is  also  included  here  in  order  to  present  a  complete  picture  of  the  variable 
medium  studies  carried  out  under  this  Contract. 
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in  which  cate  the  deflection  angle  ha*  a  maximum.  In  the  latter  case,  two  raya  con¬ 
tribute  in  any  given  direction  and  one  haa  to  diatinguiah  between  the  two  contributiona 
by  making  uae  of  the  interference  pattern.  Once  the  acattering  data  haa  been  trana- 
latad  into  a  relationahip  between  acattering  angle  and  impact  parameter,  one  may' 
proceed  to  calculate  the  electron  denaity  diatribution  in  two  di^erent  waya. 

We  can  calculate  the  deflection  angle  aa  a  function  of  the  impact  parameter 
for  an  arbitrary  radial  variation  of  electron  denaity  and  obtain  an  integral  expreaaion 
involving  the  electron  denaity  diatribution  n  which  may  be  regarded  aa  an  integral 
equation  for  n.  Thia  equation  ia  of  the  Abel  type  and  can  be  aolved.  The  final  reault 
ylelda  a  pair  of  expreaaiona  for  the  electron  denaity  and  radiua  in  terma  of  a  parameter; 
theae  conatitute  a  relationahip  between  denaity  and  radiua  in  parametric  form.  Since 
theae  e3q>reaaiona  involve  the  evaluation  of  integrala  of  producta  of  acattering  angle 
(available  in  numerical  form  from  experimental  data)  and  a  aingular  kernel,  they  are 
difficult  to  compute  numerically.  A  claaa  of  analytic  acattering  angle  functiona  waa 
introduced,  for  which  the  calculation  can  be  carried  out.  One  can  then  chooae  from 
thia  claaa  a  function  that  fita  the  experimental  data  beat  and  uae  the  correaponding 
analytical  reault  for  the  electron  denaity  diatribution. 

Alternatively,  one  can  avoid  the  aolution  of  the  integral  equation  altogether 
and  aimply  calculate  the  acattering  angle  function  for  aome  modela  of  electron  denaity 
diatribution.  One  then  chooaea  that  model  for  which  the  calculated  acatterirtg  angle 
function  fita  the  experimental  data  beat.  Three  modela  were  analyzed  in  thia  manner. 

In  all  the  caaea  of  an  underdenae  (penetrable)  column  which  were  conaidered, 
a  relationahip  between  the  maximum  acattering  angle  0^  and  the  ratio  of  maximum 
electron  denaity  in  the  column  to  the  denaity  at  which  the  plaama  becomea  impenetrable, 
waa  examined.  It  waa  found  that  in  each  caae  thia  relationahip  could  be  reaaonably 
well  approximated  by 

\  (1) 

ao  that  maximum  plaama  denaity  could  be  deduced  from  maximum  acattering  angle. 

The  detaila  of  the  above  calculationa  were  iaaued  in  a  report^  and  aubae- 
2 

quently  publiahed  . 

1.  J.  Shmoya,  *A  Propoaed  Oiagnoatic  Method  for  Cylindrical  Plaamaa*,  PIBMRI- 
828-60,  Microw.  Rea.  Inat. ,  Polytech.  Inat.  of  Brooklyn,  May  I960. 

2.  3.  Shmoya,  *A  Propoaed  Diagnoatic  Method  for  Cylindrical  Plaamaa*.  J.  Aonl. 
Phya.,  32,  pp.  689-695,  April  1961. 
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(b)  (c) 

Fig.  16  '  Ray  paths  in  cylindrically  stratified  medium 
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c.  Reflection  and  tr>inmi««ion  of  electromagnetic  wave»  in  the  preaence  of 
electron  dentity  gradients 
J.  Shmoys 

Albini  and  Jabn^  published  a  calculation  of  the  complex  reflection  and  trans¬ 
mission  coefficients  for  electromagnetic  waves  incident  normally  on  a  plane  stratified 
transition  between  free  space  and  a  lossy  isotropic  plasma.  The  variation  of  electron 
density  j^dielectric  constant  j  in  the  transition  region  was  assumed  linear  or  piecewise 
linear:  the  variation  of  magnitude  and  phase  of  the  transmission  and  reflection  coef¬ 
ficients  was  calculated  as  a  function  of  thickness  of  the  transition  layer.  Albini  and 
Jahn  concluded  that  the  coefficients  depend  strongly  on  the  width  of  the  transition  zone 
and  less'  strongly  on  the  detailed  profile  of  the  transition;  they  observed  periodicity 
in  the  amplitude  and  phase  of  the  reflection  coefficient. 

In  order  to  check  the  above  conclusions,  the  variation  of  the  same  quantities 

with  effective  layer  thickness  was  calculated  for  the  Epstein  (inverse  hyperbolic  tan- 
2 

gent)  transition  .  On  the  basis  of  the  results  obtained,  it  can  be  said  that  Albini  and 
Jahn's  conclusions  are  not  warranted.  The  magnitude  of  the  reflection  coefficient 
tends  to  zero  with  increasing  layer  thickness  much  faster  for  the  Epstein  transition, 
and  in  a  completely  different  manner.  The  phase  of  the  reflection  coefficient  does  not 
show  the  strong  periodic  component  exhibited  in  the  linear  case.  It  is  likely  that  the 
behavior  of  these  quantities  in  the  case  of  a  linear  or  piecewise  transition  is  related  to 
discontinuities  in  slope. 

It  was  originally  planned  to  issue  the  results  for  the  Epstein  transition  in  a 

memorandum  but  some  numerical  mistakes  were  found  in  the  computation  of  phase 

2 

and  it  was  felt  that  all  phase  curves  presented  in  progress  reports  should  be  recal¬ 
culated;  this  has  not  yet  been  done. 


1.  F.  A.  Albini  and  R,  G.  Jahn,  "Reflection  and  Transmission  of  Electromagnetic 
Waves  at  Electron  Density  Gradients",  J.  Appl.  Phys. ,  32,  pp.  75-82,  Jan.  1961. 

2.  10th,  11th  and  12th  Quarterly  Reports. 
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D.  Gaidlng  and  Scattering  Propertie«  of  Vnriable  Impedance  and  AbBorbing  Surfaces 
1.  Introduction 

A  reasonably  extensive  literature  exists  on  problems  of  diffraction  by 
perfectly  conducting  objects  of  simple  sLape,  describable  by  single  coordinate  surfaces 
in  an  appropriate  coordinate  system  which  renders  the  vector  or  scalar  wave  equation 
separable.  Rigorous  field  solutions  are  obtained  by  the  method  of  separation  of 
variablesi  and  their  asymptotic  evaluation  and  subsequent  physical  interpretation  in 
the  limit  of  high  and  low  frequencies  grants  considerable  insight  into  the  diffraction 
mechanism  in  these  frequency  domains.  Keller's  geometrical  theory  of  diffraction^ 
has  been  particularly  promising  in  systematizing  the  various  diffraction  processes 
which  take  place  on  relatively  arbitrary  objects  in  the  optical  (short-wavelength)  limit. 
The  validity  of  this  as  yet  non-rigorous  theory  has  been  verified  in  special  cases  by 
.:omparlng  its  predictions  with  results  obtained  from  the  above  mentioned  class  of 
exactly  solvable  problems.  The  quasi-optic  field  comprises  the  conventional  geometric- 
optical  incident  and  reflected  rays,  and  also  diffracted  rays  which  are  launched  by 
certain  special  points  on  the  object  surface  (e.  g.  t  edges,  corners,  tips,  geometrical 
shadow  boundaries,  etc.)  and  which  account  for  diffraction  effects.  Keller's  theory 
can  generally  not  predict  the  starting  amplitude  of  a  diffracted  ray  but  requires  for 
this  information  the  rigorous  solution  of  a  "canonical*  problem  which  is  locally  similar 
to  a  given  scattering  object. 

When  applying  an  asymptotic  theory  such  as  Keller's,  or  its  equivalent,  a 
point  of  importance  concerns  its  limits  of  applicability.  The  mathematical  require¬ 
ment  k-*w,  where  k  is  the  wavenumber  in  the  medium,  can  frequently  be  relaxed  and 
phrased  as  kd>  C,  where  d  is  a  characteristic  dimension  of  the  scatterer  and  C  a 
bounded  positive  constant.  For  example,  in  the  probleVn  of  diffraction  by  a  perfectly 
conducting  sphere  or  cylinder  with  radius  d,  the  first-order  asymptotic  results  are 
quite  good  when  C  10.  To  make  similar  comparisons  for  non-perfectly  conducting 
objects,  it  is  necessary  to  obtain  rigorous  solutions  for  simple  scatterers  which  ex¬ 
hibit  the  desired  surface  impedance  properties.  Seemingly  least  complicated  is  the 
specification  of  a  constant,' non-zero  surface  impedance  Z^.  It  is  found,  however, 
that  many  problems  which  are  separable  when  Z^  =  0  cease  to  be  so  when  Z  equals 
a  non-vanishing  constant.  A  well-known  example  in  this  category  is  the  perfectly 
conducting  wedge  or  half-plane;  paradoxically,  the  wedge  problem  (and  others)  separate 
for  certain  special  variations  of  Z^.  Thus,  our  quantitative  knowledge  of  diffraction 
processes  at  high  frequencies  is  much  less  complete  when  0  since  exact  solutions 
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are  scarcer  and  more  difftcult  to  obtain. 

At  the  time  of  initiation  of  this  contract,  no  rigorous  results  were  available 
for  scattering  by  a  variable  impedance  structure.  A  non-constant  surface  impedance 
is  likely  to  involve  limitations  additional  to  those  arising  from  the  object’s  structural 
shape  since  the  impedance  variation  per  local  wavelength  introduces  a  new  parameter. 
While  quasi-optic  approximations  may  be  expected  to  hold  for  slow  variations  they 
are  bound  to  be  Inaccurate  or  even  invalid  when  the  variation  is  rapid.  To  confirm 
these  expectations,  three  types  of  structures  were  investigated;  a)  a  wedge  and  a 
cone  with  a  surface  Impedance  (or  admittance,  depending  on  the  polarisation  of  the 
incident  field)  which  varies  linearly  with  distance  from  the  edge  or  tip,  respectively; 
b)  a  clrciUar  cylinder  with  a  surface  impedance  having  a  sinusoidal  peripheral  varia¬ 
tion  about  a  constant  value;  and  c)  a  plane  surface  with  a  relatively  arbitrary  (bilinear) 
variation  of  Z^.  These  choices  were  dictated  primarily  by  the  fact  that  rigorous  re¬ 
sults  are  obtainable  for  certain  two-dimensional  forms  of  excitation:  complete  solu¬ 
tions  in  category  a)  can  be  found  by  the  method  of  separation  of  variables,  while  those 
in  the  non-separable  category  b)  are  in  the  form  of  a  perturbation  series.  A  novel 
method  was  applied  to  tolve  problem  c).  The  structures  in  a)  permit  the  study  of  the 
effects  of  impedance  variation  on  plane  and  quasi-plane  surfaces,  and  near  an  edge  or 
tip,  while  those  in  b)  highlight  the  effects  of  impedance  change  on  a  curvad  object. 

Item  c)  generalises  some  of  the  results  in  a)  by  permitting  a  more  general  type  of 
Impedance  variation.  From  an  investigation  of  the  asymptotic  properties  of  these 
solutions,  carried  out  for  a)  and  b),  estimates  on  the  limits  of  validity  of  geometrical 
optics  were  obtained,  and  results  were  found  for  the  field  behavior  when  geometrical 
optics  is  no  longer  applicable.  During  the  study  of  the  problems  in  a),  a  new  type  of 
surface  wave  was  discovered  which  propagates  along  a  variable  impedance  plane,  and 
which  is  of  special  interest  in  connection  with  radiation  from  a  tapered  surface  wave 
antenna. 

To  assess  the  influence  of  a  constant  surface  impedance  on  the  field  dif¬ 
fracted  by  an  edge,  a  study  of  diffraction  by  an  imperfectly  conducting  half-plane  was 
also  carried  out,  with  emphasis  on  those  values  of  Z^  which  yield  maximum  absorption 
of  the  incident  wave  (black  screen). 

2.  Wedge  or  cone  with  a  linearly  varying  surface  impedance 

This  study  deals  with  a  line-source  excited  wedge  (Fig.  17(a)  )  or  a  ring- 
source  excited  cone  (Fig.  17(b))  whose  surface  impedance  (or  admittance)  varies  lin¬ 
early  with  distance  from  the  apex  or  tip,  respectively.  For  the  wedge  problem,  two 
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different  impedance  valuea  and  are  permitted  on  the  facet  at  ^  ■  0  and  ^  ■  a, 
reapectively;  when  the  exciting  line  source  comprises  magnetic  currentSi  one  has 
Z  B  ikpc  •  where  c  and  c  are  constants  subject  to  the  restriction  Imc.  .<0 
to  assure  passiT^ty.  k  is  the  real  wavenumber  in  the  surrotmding  medium  and  p  is  the 
radial  variable.  Analogous  considerations  apply  to  Z^  in  the  cone  problem  except 
that  p  la  replaced  by  the  spherical  variable  r.  These  boundary  value  problems  can 
be  solved  by  the  method  of  separation  of  variables)  and  alternative  representations 

which  are  especially  suited  to  the  study  of  the  various 


=0 


(a)  WEDGE  (b)  SEMI -INFINITE  CONE 

Fig.  17  -  Physical  structures 

propagation  and  diffraction  phenomena. 

a.  Formulation  err.phasizing  surface  guiding  properties 

If  the  wedge  in  Fig.l7(a)  is  viewed  as  a  radial  waveguide  in  which  cylindri- 
cally  spreading  or  contracting  waves  propagate)  the  appropriate  representation  in¬ 
volves  the  complete  set  of  (discrete)  angular  eigenfunctions.  During  the  investigation 
of  the  angular  eigenvalue  problem)  a  new  type  of  surface  wave  was  discovered^  which 
is  guided  by  one  or  both  of  the  wedge  faces,  propagates  outward  from  the  apex  like  a 
cylindrical  wave)  and  decays  in  the  ^-direction  away  from  the  wedge  surface.  Its 
energy  in  the  region  between  the  wedge  face  at  ^  >  0  and  a  plane  parallel  to  4  ■  0,  for 
example)  does  not  remain  constant;  instead,  continuous  radiation  takes  place  because 


of  the  solution  were  constructed 


<f>sa 
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of  the  surface  impedance  variation.  This  wave  is  particularly  interesting  since  it 

may  serve  as  a  rigorous  prototype  for  the  analysis  of  a  tapered  surface  wave  antenna. 

2 

A  theoretical  study  of  the  latter  structure  was  carried  out  and  it  was  shown  that  very 
high  endfire  gains  can  be  achieved  if  a  constant  reactance  surface  waveguide  is  ter¬ 
minated  by  a  suitable  tapered  reactance  section  from  which  the  radiation  takes  place. 
The  calculations  also  included  data  on  the  sidelobes  to  be  expected  from  the  junction 
between  the  two  waveguide  sections. 

The  field  representation  in  terms  of  radially  propagating  waves  is  also 
useful  for  the  study  of  the  diffraction  fields  near  the  apex.  Analogous  considerations 
apply  to  the  cone  problem. 

b.  Formulation  emphasizing  scattering  properties 

To  analyze  the  scattering  properties  of  the  wedge  (or  cone)  when  both  the 
source  and  observation  points  are  located  many  wavelengths  from  the  apex,  an  angular 
transmission  representation  is  appropriate.  From  this  viewpoint,  waves  are  regarded 
to  propagate  along  the  41  (or  0)  directions,  and  the  field  representation  is  in  terms  of 
the  (continuous)  set  of  radial  eigenfunctions.  In  its  most  elementary  form,  the  angular 
space  is  infinitely  extended,  and  the  effects  of  boundaries  delimiting  the  finite  physical 
^  (or  6)  space  are  taken  into  account  by  multiple  reflection  or  image  construction.  The 
general  utility  of  the  angular  transmission  approach  in  the  analysis  of  high-frequency 
diffraction  problems  was  discussed  in  reference  3.  A  first-order  asymptotic  evalua¬ 
tion  of  the  formal  integral  solution  yields  a  result  for  the  far  field  which  can  be  inter¬ 
preted  conveniently  in  terms  of  the  geometric-optical  ray  picture  shown  in  Fig.  18. 


Fig.  18  -  Quasi-optic  fields 
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The  field  eomprleee  direct  and  reflected  raye  -  each  in  a  domain  predicted 
from  geometrical  optica;  in  addition,  a  cylindrical  diffracted  wave  emanating  from  the 
edge  ia  obaervable  in  both  the  illuminated  and  ahadow  regiona.  The  reflected  ray 
amplitude  la  determined  by  the  local  wedge  impedance  at  the  point  of  reflection,  thereby 
confirming  the  validity  of  ita  conatruction  by  geometric -optical  conaiderationa.  The 
diffracted  wave  amplitude  D  la  found,  however,  to  depend  not  only  on  the  impedance  at 
.the  diffraction  point  (the  edge)  but  ia  influenced  by  ita  rate  of  variation  in  the  vicinity 
of  the  edge.  If  the  aurface  impedance  variation  at  ^  -  0  ia  written  aa  Ap,  and  the  wedge 
face  at  ^  «  a  ia  aaaumed  to  be  perfectly  conducting,  then  when  A«  1,  D  ia  eaaentially 
the  aame  aa  for  a  zero  impedance  wedge;  however,  when  A»  1,  O  ia  eaaentially  the 
aame  aa  for  a  wedge  with  infinite  impedance  at  s  0  although  the  im)pedance  at^the  edge 
p  «  0  ia  zero,  A«1  and  A»  1  correapond,  reapectively,  to  alow  and  rapid  impedance 
varlationa  near  the  edge  whence  it  ia  noted  that  the  geometric -optical  predictiona  apply 
approximately  for  alow,  but  are  invalid  for  rapid,  variation.  The  behavior  of  D  for 
intermediate  rangea  of  A  haa  been  calctilated  in  reference  2  which  alao  containa  a  com¬ 
plete  aaymptotic  expanaion  of  the  diffraction  field. 

3.  Circular  cylinder  with  peripherally  varying  surface  impedance 
C.  J.  Marcinkowski  and  L.  B.  Felsen 


Thia  aection  deala  with  the  diffraction  of  a  line  aource  current  field  by  a 
cylinder  with  the  ainuaoidally  varying  aurface  impedance  ^ 


1  +  a  coa 


P('t»-'l>o)] 


(1) 


where  ^  ia  the  angle  of  obaervation,  ia  an  arbitrary  phaae  angle,  and  p  expreaaea 
the  periodicity  of  the  aurface  impedance.  Thia  atudy  waa  undertaken  for  the  purpoae 
of  exploring  the  aimultaneoua  influence  of  aurface  curvature  and  variability  of  aurface 
propertiea  on  the  field  diffracted  by  an  object  at  high  frequenciea.  The  ainuaoidal 
impedance  variation  ia  highly  deairable  aince  an  appropriate  choice  of  the  periodicity 
parameter  p  permita  the  aimulation  of  alow,  monotonic  varlationa  on  the  one  hand, 
and  of  rapid  fluctuationa  on  the  other.  For  j a|  « 1,  the  problem  can  be  aolved  by  a 
perturbation  method  which  leada  to  a  Green'a  function  repreaentation  patterned  on  the 
well  known  reaulta  for  a  conatant  impedance  cylinder.  The  infinite  angular  trana- 
miaalon  line  approach  mentioned  in  Sec.  2b  waa  alao  employed  here  to  derive  aolutlona 
for  the  conatant  impedance  cylinder  which  permit  a  direct  aaymptotic  field  evaluation 
in  both  the  illuminated  and  ahadow  regiona.  Theae  reaulta  were  then  extended  to  in¬ 
clude  the  variable  impedance  caae  for  which  the  fielda  were  expreaaed  aa  power  aeriea 
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involving  the  perturbation  parameter  a. 

These  formal  solutions  were  then  investigated  to  first  order  in  a  for  a 
plane  wave  incident  on  a  large  cylinder.  ^  A  saddle  point  evaluation  yielded  a  result 
which  is  interpretable  in  terms  of  the  geometrical  optics  appropriate  to  a  cyllndrl- 
cally  curved,  convex  reflection  grating.  The  solution  comprises  the  first  and  higher 
order  reflected  rays  familiar  from  the  theory  of  the  plane  grating,  modified  by  ap¬ 
propriate  divergence  coefficients.  ^  This  extension  of  conventional  geometrical  optics 
provides  a  complete  and  qviantltatlve  physical  picture  of  the  mechanism  of  reflection 
of  waves  from  a  curved  surface  with  periodic  Impedance  variation. 

For  a  slowly  varying  surface  Impedance  (p  « 1),  conditions  were  derived 
under  which  the  various  reflected  rays  in  the  solution  to  order  a  may  be  combined 
into  a  single,  specularly  reflected  ray  with  a  reflection  coefficient  depending  only  on 
the  local  impedance  value  at  the  point  of  reflection.  ^  Since  this  is  precisely  the  be¬ 
havior  of  the  geometric -optical  field  reflected  from  a  variable,  non-periodic  impedance 
surface  these  conditions  which  involve  both  the  rate  of  change  of  surface  impedance  and 
the  radius  of  curvature  of  the  surface,  serve  as  criteria  for  the  validity  of  geometrical 
optics  when  applied  to  curved,  variable  impedance  structures. 

4.  Plane  surface  with  bilinear  impedance  variation 

C.  J.  Marcinkowski 

The  problem  of  diffraction  by  a  linearly  varying  surface  impedance  on  a 
wedge  and  cone  was  solved  previously  by  the  method  of  separation  of  variables  (see 
Sec.  2).  To  deal  with  the  problem  of  a  linearly  varying  impedance  on  an  infinite  plane 
surface,  separation  of  variables  is  Inapplicable  and  other  methods  of  analysis  must  be 
employed.  The  Impedance  function  is  assumed  to  have  the  form 

2:{x)  =  y^{to  +  .  -'»<x<«»  ,  (2) 

where  k  is  the  free-space  wavenumber,  the  free-space  impedance,  and  and 

are  arbitrary  complex  numbers.  By  a  combination  of  Fourier  integral  and  integra- 
tion-by-parts  procedures,  the  problem  for  the  unknown  Fourier  transform  is  reduced 
to  a  first-order  differential  equation  which  was  solved  explicitly  by  a  standard  proce¬ 
dure.^  Since  this  method  of  solution  has  apparently  not  been  considered  before,  an 
ei^loratlon  of  this  technique  was  carried  out  to  determine  the  extent  to  which  It  could 
be  developed.  This  led  to  the  discovery  of  the  formal  solutions  for  a  surface  Impedance 
in  the  bilinear  form^ 
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Z(x) 


to  + 

-f  a,  kx 

O  X 


where  and  are  arbitrary  complex  numbers. 


(3) 


The  solution  for  the  linear  impedance  (2)  in  the  neighborhood  |  « 1  is 

of  physical  importance  since  it  involves  the  interesting  quasi-optical  region  of  a 

slowly  varying  surface  impedance.  ^  This  neighborhood  is  also  of  mathematical 

Importance  because  the  solution  contains  the  factor  l/^,  in  an  exponential  and  as  a 

10  ^ 

multiplier  in  front  of  the  exponential,  thereby  exhibiting  an  apparent  singularity  at 
«  0.  For  these  reasons  the  limiting  behavior  |  C J  -*  0  was  investigated  for  the 

i  *  7 

physically  important  case  of  a  reactive  surface  impedance.  It  was  found  that  solu¬ 
tions  were  allowed  for  a  constant  inductive  reactance  arg((^)  ~ --ir/2,  and  for  either 
type  of  varying  reactance  with  arg((p  at  ±  v  fZ.  Dual  results  were  obtained  for  surface 
admittance  functions.  This  produced  the  formal  solutions  for  the  infinite  plane  with 
the  blllnearly  varying  surface  impedance^  but  no  time  was  available  for  a  detailed 
examination  of  the  physical  properties  of  any  of  these  solutions. 


5.  Diffraction  by  an  absorbing  half-plane 
C.  J.  Marcinkowski 


Although  the  two-dimensional  electromagnetic  problem  of  diffraction  by  a 
constant  impedance  half-plane  has  been  studied  by  previous  workers,  there  still  re¬ 
main  two  major  deficiencies  in  the  studies  reported  in  the  literature.  One  of  these 
lies  in  the  mathematical  complexity  of  the  formal  solution  while  the  other  is  due  to 
the  very  narrow  range  of  surface  impedance  values  which  have  been  examined.  Several 
steps  were  taken  to  simplify  the  solution.  The  first  was  to  make  use  of  a  much  simpler 
but  hitherto  neglected  result  obtained  by  Grinberg  and  Fock  and  used  extensively  by 
the  Russian  school  of  investigators.  The  surface  impedance  was  described  in  terms 
of  an  angle  of  "perfect  absorption"  designated  as  ,  the  angle  at  which  an  incident 
plane  wave  produces  no  reflected  plane  wave.  By  means  of  this  procedure,  the  asymp¬ 
totic  far-field  scattering  patterns  could  be  specified  in  a  physically  significant  manner 

g 

in  terms  of  three  similar,  dimensionless  parameters  ;  the  angle  of  incidence  ,  the 
angle  of  observation  and  the  angle  of  perfect  absorption 

Since  previous  investigations  have  been  restricted  to  the  narrow  range  of 
snmll  surface  impedances,  they  have  supplied  little  information  concerning  the  in¬ 
fluence  of  surface  impedance  on  the  scattered  far  fields.  In  contrast,  the  scattered 
far  fields  were  studied  under  this.  Contract  for  a  relatively  arbitrary,  constant  surface 
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impedance  which  may  act  like  a  black  screen.  A  null  was  discovered  in  the  back- 
scattered  fields  with  properties  which  permitted  the  interpretation  of  the  asymmetric 
patterns  of  the  constant  impedance  half-plane  in  terms  of  an  angular  deformation  or 

g 

distortion  of  the  more  symmetric  patterns  of  a  perfectly  reflecting  structure.  The 
scattering  patterns  were  discovered  to  be  roughly  similar  to  those  of  the  Sommerfeld 
black  screen®  which  has  the  very  simple  amplitude  behavior  ^  +  +£  "  • 

All  these  simplifications  permit  us  to  sketch  oat  the  tnain  properties  of  the  diffracted 

g 

far  fields  in  a  simple  qualitative  manner  without  the  need  for  elaborate  calculations. 
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E.  Other  Pllfraction  Problema 

1.  Plane  wave  back-»cattering  from  a  temi-infinite  cone 

In  previous  studiee  of  diffraction  by  a  eemi -infinite  cone\  it  haa  been 

shown  that  the  far  field  can  be  decomposed  into  incident,  reflected  and  diffracted 

constituents.  The  first  two  are  the  geometric -optical  field  contributions  which  can 

be  evaluated  in  closed  form  while  the  third  involves  a  canonical  diffraction  integral 

which  specifies  the  angular  variation  of  the  spherical  wave  emanating  from  the  cone 

tip.  Since  an  exact  evaluation  of  the  diffraction  integral  does  not  seem  possible,  it 

was  decided  to  carry  out  an  approximate  evaluation  for  the  special  case  of  plane  wave 

back-scattering  along  the  cone  axis.  Approximate  results  had  previous.’.y  been  ob- 

2 

tained  for  the  limiting  cases  of  cones  having  very  small  or  large  apex  angles  ;  the 
present  study  yielded  a  formula  applicable  for  arbitrary  angles,  and  the  resulting 
magnitude  of  the  back- scattered  field  has  been  given  in  the  1st  Quarterly  Report  (see 
also  reference  3). 

2.  Relation  between  a  class  of  two-dimensional  and  three-dimensional  diffraction 
problems 

By  means  of  a  certain  transformation,  a  relationship  was  demonstrated 
between  a  class  of  two-dimensional  and  three-dimensional  scalar  or  electromagnetic 
diffraction  problems.  The  basic  three-dimensional  configuration  consists  of  a  per¬ 
fectly  reflecting  half-plane  excited  by  a  ring  source  centered  about  the  edge  and 
having  a  variation  exp(J;  i^/2),  where  is  the  azimuthal  variable;  in  addition,  a 
perfectly  reflecting  rotationally  symmetric  obstacle  whose  surface  is  defined  by 
f(p,  z)  s  0  (p,  z  are  cylindrical  coordinates),  may  be  superposed  about  the  edge  (z-axis). 
This  problem  was  shown  to  be  simply  related  to  the  two-dimensional  one  for  the  line 
source  excited  configuration  f(y,  z)  a  0,  where  y  and  z  are  Cartesian  coordinates  (see 
^8<19)*  Various  special  obstacle  configurations  were  treated  in  detail. 

For  the  general  case  of  arbitrary  electromagnetic  excitation,  the  above- 
mentioned  transformation  was  used  to  construct  the  solution  for  the  diffraction  by  a 
perfectly  conducting  half-plane  from  the  knowledge  of  appropriate  scalar  solutions. 
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Fig.  19  -  Related  configurations 


namely  those  which  obey  the  same  equations  and  boundary  conditions,  and  have  the 
same  excitations,  as  the  Cartesian  components  of  the  electromagnetic  field. 

The  results  of  this  analysis  have  been  reported  in  reference  1. 

F.  Cerenkov  Radiation 

L.  B.  Felsen  and  A.  Hessel 

Cerenkov  radiation  problems,  dealing  with  radiation  from  a  charge  (or 
I  charges)  in  straight  uniform  motion  near  or  in  perturbing  media  which  can  support 
electromagnetic  waves  with  phase  speeds  slower  than  the  particle  speed,  have  been 
analyzed  in  the  literature  by  various  methods  which  become  quite  involved  when  the 
fields  are  not  derivable  from  a  single  scalar  potential  function.  For  example,  a  charge 
moving  parallel  to  the  plane  boundary  of  a  plane  stratified  dielectric  medium,  a  plane 
grating,  a  plane  reactive  surface,  etc. ,  excites  both  E  and  H  modes  in  the  z-direction 
perpendicular  to  the  interface  and  the  associated  field  problem  is  a  vector  problem. 
While  modal  network  procedures  have  been  applied  for  some  time,  and  with  considerable 
success,  to  the  analysis  of  radiation  from  stationary  sources,  they  had  not  been  utilized 
for  fields  excited  by  moving  sources.  Such  a  study  was  irdtiated  under  this  contract 

1.  L.  B.  Felsen  and  S.  N.  Karp,  "Relation  between  a  Class  of  Two-Dimensional  and 
Three-Dimensional  Diitfraction  Problems",  Report  R-694-58,  PIB-622,  Microw. 

Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Jan.  1959. 
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(•ee  9th  Quarterly  Report)*  and  it  waa  ahown  how  the  aolution  of  Cerenkov  radiation 
problema  can  be  ayatematized  through  the  uae  of  modal  techniquea  by  which  one  decom- 
poaea  the  temporal  Fourier  tranaform  of  the  particle  field  into  E  and  H  modea  along  z. 
The  z-dependent  modal  amplitudea  are  evaluated  from  a  network  problem  wherein  the 
perturbing  interface  or  structure  is  represented  by  an  equivalent  network  which  gener¬ 
ally  couples  the  various  modes,  and  the  time  dependent  fields  are  recovered  by 
Fourier  synthesis. 


For  example,  if  a  charge  moves  parallel  to  a  periodic  structure  (e.  g. ,  a 
grating)  (see  Fig. 20(a)  ),  the  associated  network  problem  is  the  one  schematized  in 
Fig.  20(b).  The  temporal  Fourier  transform  of  the  current  excitation  descriptive  of 
the  moving  charge  is  a  line  source  with  a  linearly  progressing  phase.  This  source 
drives  the  (n  =  0,  non-propagating)  £  and  H  mode  transmission  lines,  and  coupling  to 
the  propagating  mode  lines  denoted  by  subscripts  n  to  ni- N  is  achieved  by  the  coupling 
network  representative  of  the  grating.  Since  the  source  moves  in  vacuum  where  the 
particle  speed  is  less  than  the  velocity  of  light,  the  associated  fields  are  exponentially 
decaying  and  radiation  occurs  only  because  of  the  coupling  to  one  or  more  of  the  propa¬ 
gating  grating  modes.  . 
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(b)  EQUIVALENT  NETWORK 


Fig.  20  -  Cerenkov  radiation 
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The  above-described  method  has  been  presented  in  reference  1,  and  various 
applications  have  been  under  study.  The  results  will  be  reported  under  the  successor 
contract. 

G.  Network  Theory 

D.  C.  Youla 

The  amount  of  material  reported  in  this  period  is  varied  and  quite  extensive. 

In  the  first  quarterly,  July  1958  -  September  1958,  appeared  two  notes 
dealing  with  parametrically  excited  structures.  The  first,  ‘Response  of  Parametri¬ 
cally  Excited  Networks",  by  L.  Smilen  and  D.  C.  Youla,  developed  a  single-frequency 
network  formalism  enabling  the  analysis  of  cascades  of  passive  networks  with  em¬ 
bedded  varactors  to  be  carried  through  in  relatively  simple  fashion.  The  second, 

*A  Generalization  of  the  Manley-Rowe  Relations",  by  D.  C.  Youla,  developed  power 
frequency  formulas  for  non-hysteretic  lossy  varactors.  The  main  conclusion  was  the 
fact  that  formulas,  unlike  the  lossless  case,  depended  strongly  on  the  varactor  charac¬ 
teristic. 

In  the  second  quarterly,  October  1958  -  December  1958,  a  note  by  D.  C.  Youla 
entitled  "Synthesis  of  Smooth  Non-Uniform  Transmission  lines  from  a  Prescribed 
Scattering  Coefficient",  presented  an  exact  technique  for  the  synthesis  of  a  smooth, 
lossless  non-uniform  transmission  line  (of  finite  length)  from  prescribed  input  reflec¬ 
tion  datum.  The  method  reduced  the  problem  to  the  solution  of  a  one-parameter 
family  of  Fredholm  integral  equations.  As  yet,  a  truly  significant  engineering  adap¬ 
tation  is  not  available  and  a  good  deal  remains  to  be  done  along  these  lines. 

In  the  third  quarterly,  January  1959  -  March  1959,  the  article,  "A  Note  on 
the  Stability  of  Linear,  Non-Reciprocal  N-Ports",  by  D.  C.  Youla,  suggested  a  method 
for  determining  the  real-frequency  stability  of  a  non-reciprocal  linear,  time-invariant 
n-port  via  its  impedance  matrix.  Inherent  limitations  as  well  as  advantages  were 
pointed  out.  An  application  to  the  2-port  case  immediately  yielded  the  classical 
formulas  due  to  Ralsbeck. 

The  fourth  quarterly,  April  1959  -  June  1959,  contains,  in  the  paper, 
"Weissfloch  Equivalents  for  Lossless  Zn-Ports"  by  O.  C.  Youla,  a  rather  significant 
contribution  to  the  representation  of  lossless  multi-mode  discontinuities.  Work  on 

1.  L.  B.  Felsen  and  A.  Hessel,  "A  Network  Approach  to  the  Analysis  of  Cerenkov 
Radiation  Problems',  Electrophysics  Memorandum  No.  60,  PIBMRI-878-60, 

Microw.  Res.  Inst. ,  Polytech.  Inst,  of  Brooklyn,  Nov.  9,  I960. 
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associated  measurement  techniques  is  now  in  progress. 

In  'Synthesis  of  Networks  Containing  both  Positive  and  Negative  Resistors" > 
by  D.  C.  Youlai  fifth  quarterly,  July  1959  -  September  1959,  a  network  theory  of  the 
ustial  lumped  elements  plus  negative  resistors  is  initiated.  It  is  shown  that  many  of 
the  passive  network  techniques  and  theorems  may  be  used  to  derive  results  for  the 
active  case.  Two  specific  theorems  concerning  the  number  of  positive  and  negative 
resistors  required  for  the  synthesis  of  lumped  active  n-ports  are  derived  and  their 
implications  discussed  in  some  detail. 

In  the  sixth  quarterly,  October  1959  -  December  1959,  a  rather  complete 
analysis  was  presented  of  the  single-frequency  invariantive  properties  of  a  linear, 
time -invariant  n-port.  In  the  note,  'The  Analytic  Foundations  of  the  "Spot"  Frequency 
Theory  of  Linear,  Noiseless,  Time -Invariant  n-Ports",  by  D.  C.  Youla,  the  concept 
of  the  cross -ratio  matrix  of  four  n-ports  was  introduced  for  the  first  time.  This 
notion  led  naturally  to  the  'self"  cross-ratio  matrix  of  a  single  non-reciprocal  n-port. 
A  most  important  conclusion  was  that  the  ei^genvalues  of  this  'self"  cross-ratio  ma¬ 
trix  remained  invariant  under  lossless,  reciprocal  2n-port  embedding.  That  these 
eigenvalues  reflected  everything  of  physical  interest  was  demonstrated  by  exhibiting 
several  canonical  forms. 

The  note,  'The  Inverse  Problem  for  Smooth,  Lossless  Non-Uniform  Trans¬ 
mission  Lines',  by  D.  C.  Youla,  seventh  quarterly,  January  1960  -  March  I960,  con¬ 
tinues  the  related  research  reported  in  the  second  quarterly  by  the  same  author. 
Additional  theorems  and  refinements  are  supplied  but  the  goal  of  a  practical,  exact 
engineering  synthesis  technique  is  not  attained  and  further  work  is  necessary. 

Many  problems  in  electrical  engineering,  such  as  the  synthesis'^of  linear 
n-ports  and  the  detection  and  filtration  of  multivariable  systems  corrupted  by  sta¬ 
tionary  additive  noise  depend  for  their  successful  solution  upon  the  factorization  of  a 
matrix-valued  function  of  a  complex  variable  p.  The  paper,  "On  the  Factorization  of 
Rational  Matrices",  by  D.  C.  Youla,  eighth  quarterly,  April  I960  -  June  I960,  presents 
several  algorithms  for  affecting  such  decompositions  for  the  class  of  rational  matrices 
G(p);  i.  e. ,  matrices  whose  entries  are  ratios  of  polynomials  inp.  The  methods  em¬ 
ployed  are  elementary  in  nature  and  center  around  the  Smith  canonical  form  of  a 
polynomial  matrix.  Several  nontrivial  examples  are  worked  out  in  detail  to  illustrate 
the  theory. 
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In  "Some  Reeulta  in  the  Single-Frequency  Theory  of  Linear  Noieelesa  Net- 
worka*.  by  O.  C.  Youla,  ninth  quarterly,  July  I960  -  September  I960,  the  work  ini¬ 
tiated  in  the  aixth  quarterly  ia  generalized  and  enlarged.  The  *Q"  matrix  of  an  n-port 
ia  introduced  and  ahown  to  alwaya  exiat  irreapective  of  the  activity  character  of  the 
atructure.  The  requirementa  of  reciprocity,  paaaivity,  loaaleaaneaa  and  aymmetry 
are  formulated  in  terma  of  Q.  All  croaa-ratio  matricea  are  expreaaed  in  terma  of 
Q  which  emergea  as  an  invariantive  repreaentation  of  the  n-port.  Its  tranaformation 
propertiea  under  2n-port  embedding  are  remarkably  simple  and  completely  avoid  the 
matrix  bilinear  operators  characteristic  of  the  immittance  and  scattering  bases. 

In  the  tenth  quarterly,  October  I960  -  December  I960,  a  major  breakthrough 
in  2-port  insertion-loss  filter  design  is  reported  in  "Darlington  Synthesis  via  Richard's 
Theorem",  by  O.  C.  Youla.  Explicit  formulas  for  the  element  values  of  the  types  A, 

B,  C  and  D  sections  are  derived  and  summarized  in  chart  form  in  "A  New  Theory  of 
Cascade  Synthesis",  by  D.  C.  Youla,  Research  Report  No.  PIBMRI-916-61,  May  29, 
1961.  By  means  of  the  chart,  the  design  of  a  filter  may  be  accomplished  with  a  mini¬ 
mum  of  numerical  labor.  In  addition,  a  new  theorem  on  positive-real  functions  is 
presented  which  serves  as  the  basis  for  a  new  technique  of  1-port  cascade  synthesis. 
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Report!  and  memoranda  published  under  the  contract: 

E.  Arbel,  'Radiation  from  a  Point  Source  in  an  Anisotropic  Medium*,  PIBMRI-86I-6O1 
Microw.  Res.  Inst.,  Polytech.  Ipst.  of  Brooklyn,  Nov.  2,  I960 

E.  Arbel  and  L.  B.  Felsen,  *On  Electromagnetic  Green's  Functions  for  Uniaxially 
Anisotropic  Regions",  PIBMRI-985>61,  Microw.  Res.  Inst.,  Polytech. 

Inst,  of  Brooklyn,  Feb.  26,  1962 

E.  Arbel  and  L.  B.  Felsen,  'On  the  Theory  of  Radiation  from  Sources  in  Plane 

Stratified  Anisotropic  Plasma  Media",  PIBMRI-1069-62,  Microw.  Res. 
Inst.,  Polytech.  Inst,  of  Brooklyn,  Sept.  25,  1962 

L.  B.  Felsen  and  S.  N.  Karp,  'Relation  Between  a  Class  of  Two-Dimensional  and 

Three-Dimensional  Diffraction  Problems",  Report  R-694-58,  PIB  622, 
Microw.  Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Jan.  12,  1959 

L.  B.  Felsen,  "On  the  Electromagnetic  Properties  of  Wedges  and  Cones  with  a  Linearly 
Varying  Surface  Impedance*,  Report  R-736-59,  PIB  664,  Microw.  Res. 
Inst.,  Polytech.  Inst,  of  Brooklyn,  April  14,  1960 

L.  B.  Felsen,  'Diffraction  by  Objects  in  a  Certain  Variable  Plasma  Medium',  Report 
R-765-59,  PIB  693,  (Electrophysics  Memo.  No.  55),  Microw.  Res.  Inst., 
Polytech.  Inst,  of  Brooklyn,  July  20,  1959 

L.  B.  Felsen  and  A.  Hessel,  "A  Network  Approach  to  the  Analysis  of  Cerenkov  Radia¬ 
tion  Problems',  PIBMRI-878-60,  (Electrophysics  Memo.  No.  60),  Microw. 
Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Nov.  9,  I960 

L.  B.  Felsen,  'Radiation  from  a  Uniaxially  Anisotropic  Plasma  Half-Space",  PIBMRI- 
1058-62,  Microw.  Res.  Inst. ,  Polytech.  Inst,  of  Brooklyn,  Aug.  27,  1962 

L.  B.  Felsen,  "On  Diffraction  by  Objects  in  a  Uniaxially  Anisotropic  Medium", 

PIBMRI-1073-62,  (Electrophysics  Memo.  No.  83),  Microw.  Res.  Inst., 
Polytech.  Inst,  of  Brooklyn,  Aug.  31,  1962 

A.  Hessel,  N.  Marcuvitz  and  J.  Shmoys,  "Scattering  and  Guided  Waves  at  an  Inter¬ 
face  Between  Air  and  a  Compressible  Plasma",  PIBMRI-953-61,  Microw. 
Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Sept.  29,  1961 

C.  J.  Marcinkowski,  "Diffraction  by  an  Absorbing  Half-Plane",  Report  R-750-59, 

PIB  678,  Microw.  Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Sept.  29,  1959 

C.  J.  Marcinkowski  and  L.  B.  Felsen,  "Diffraction  by  a  Cylinder  with  a  Variable 
Surface  Impedance",  PIBMRI-821-60,  Microw.  Res.  Inst. ,  Polytech, 

Inst,  of  Brooklyn,  Jan.  27,  1961 

A.  A.  Oliner  and  T.  Tamir,  "Backward  Waves  on  Isotropic  Plasma  Slabs",  PIBMRI- 
930-61,  (Electrophysics  Memo.  No.  66),  Microw.  Res.  Inst.,  Polytech. 
Inst,  of  Brooklyn,  July  5,  1961 

A.  A.  Oliner  and  T.  Tamir,  "Radiation  from  Semi-Infinite  Slot-Excited  Plasma  Sheath 
Configurations",  PIBMRI-1029-62,  Microw.  Res.  Inst.,  Polytech.  Inst, 
of  Brooklyn,  May  17,  1962 
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Report*  and  memoranda  (continued) 

J.  Shmoye,  "A  Propoeed  Oiagnoatic  Method  for  Cylindrical  Plaamae''«  PIBMRI-828-60, 
Microw.  Rea.  Inat.  i  Polytech.  Inat.  of  Brooklyn,  May  16,  I960 

J.  Shmoya  and  P.  Hirach,  "Surface  Wave*  on  Plaama  Slab**,  PIBMRI-1093x62, 

(Electrophyaica  Memo.  No.  85),  Microw.  Rea.  Inat.,  Polytech.  Rea.  Inat., 
Oct.  18.  1962 


T.  Tamir,  "The  Electromagnetic  Field  of  a  Source-Excited  Plaama  Slab",  PIBMRI- 
932-61,  Microw.  Rea,  Inst.,  Polytech.  Inst,  of  Brooklyn,  Sept.  1961 

T.  Tamir  and  A.  A.  Oliner,  "The  Influence  of  Complex  Waves  on  the  Radiation  Field 
of  a  Slot-Excited  Plasma  Layer",  PIBMRI-959-61,  Microw.  Res.  Inat., 
Polytech.  Inst,  of  Brooklyn,  Nov.  2,  1961 

T.  Tamir  and  A.  A.  Oliner,  "The  Spectrum  of  Electromagnetic  Wayes  Guided  by  a 

Plasma  Layer",  PIBMRI-970-61,  Microw.  Res.  Inst. ,  Polytech.  Inat.  of 
Brooklyn,  Dec.  13,  1961 

O.  Youla,  "A  Generalization  of  the  Manley-Rowe  Relations",  R-702-58,  PIB  631, 

(Networks  and  Waveguide  Memo.  No.  16),  Microw.  Res.  Inat. ,  Polytech. 
Inat.  of  Brooklyn,  Dec.  17,  1958 

D.  Youla,  "Single  Frequency  Synthesis  from  a  Prescribed  Scattering  Matrix",  R-721-59, 
PIB  649,  (Networks  and  Waveguide  Memo.  No.  20),  Microw.  Res.  Inst., 
Polytech.  Inst,  of  Brooklyn,  March  15,  1959 

D.  Youla,  "Weissfloch  Equivalents  for  Lossless  2  N-Ports",  R-753-59,  PIB  681, 

(Networks  and  Waveguide  Memo.  No.  25),  Microw.  Res.  Inst. ,  Polytech. 
Inst,  of  Brooklyn,  Sept.  4,  1959 

D.  Youla,  "Some  New  and  Useful  Matrix  Results',  PIBMRI-807-60,  (Networks  and 
Waveguide  Memo.  No.  34),  Microw.  Res.  Inst,,  Polytech.  Inst,  of 
Brooklyn,  March  7,  1960 

D.  Youla,  "On  the  Factorization  of  Rational  Matrices",  PIBMRI-855-60,  Microw. 

Rea.  Inst. ,  Polytech.  Inst,  of  Brooklyn,  Sept.  20,  I960 

D.  Youla,  "A  Normal  Form  for  a  Matrix  under  the  Unitary  Congruence  Group", 

PIBMRI-858-60,  (Networks  and  Waveguide  Memo.  No.  41),  Microw.  Res. 
Inat. ,  Polytech.  Inst,  of  Brooklyn,  Aug.  30,  I960 

D.  Youla,  "Darlington  Synthesis  via  Richards  Theorem",  PIBMRI-901-61,  (Networks 
and  Waveguide  Memo.  No.  50),  Microw.  Res.  Inst. ,  Polytech.  Inst,  of 
Brooklyn,  March  14,  1961 

O.  Youla,  "A  New  Theory  of  Cascade  Synthesis",  PIBMRI-916-61,  Microw.  Res.  Inst., 
Polytech.  Inst,  of  Brooklyn,  May  29,  1961 

D.  Youla,  "Cascade  Synthesis",  PIBMRI-992-62,  (Networks  and  Waveguide  Memo. 

No.  62),  Microw.  Res.  Inst.,  Polytech.  Inst,  of  Brooklyn,  Feb.  14,  1962 

O.  Youla,  L.  Kaplan  and  D.  Stock,  "Analytical  Foundations  of  Linear  Time-Invariant 
N-Ports",  PIBMRI-1028-62,  Microw.  Res.  Inst.,  Polytech.  Inst,  of 
Brooklyn,  May  3,  1962 
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Reporf  and  memoranda  ( c ontinued) 

D.  Youla,  *  Limits  of  Errors  in  Mismatched  Attenuators,  Parti",  PIBMRI-1042-62, 
(Networks  and  Waveguide  Memo.  No.  68),  Microw.  Res.  Inst. ,  Polytech. 
Inst,  of  Brooklyn,  June  19,  1982 


Papers  published  in  scientific  journals: 

E.  Arbel  and  L.  B.  Felsen,  "On  the  Theory  of  Radiation  from  Sources  in  Anisotropic 
Plasmas,  Parts  I  and  II",  to  be  published  in  the  Proceedings  of  the  1962 
Copenhagen  Symposium  on  Electromagnetic  Theory  and  Antennas 

L.  B.  Felsen,  "Electromagnetic  Properties  of  Wedge  and  Cone  Surfaces  with  a 

Linearly  Varying  Surface  Impedance",  IRE  Trans.,  AP-7,  231  (1959) 

L.  B.  Felsen  and  S.  N.  Karp,  "Relation  Between  a  Class  of  Two-Qimensional  and 

Three-Dimensional  Diffraction  Problems",  IRE  Trans.,  AP-8,  407  (I960) 

L.  B.  Felsen,  "Radiation  from  a  Tapered  Surface  Wave  Antenna",  IRE  Trans.  PGAP, 
(I960) 

L.  B.  Felsen  and  A.  Hessel,  "A  Network  Approach  to  the  Analysis  of  Cerenkov 
Radiation  Problems",  Nuovo  Cim.  -  Serie  X,  19,  1065,  (1961) 

L.  B.  Felsen,  "High-Frequency  Diffraction  by  a  Wedge  with  Linearly  Varying  Surface 
Impedance",  Appl.  Sci.  Res.,  B9,  170  (1961) 

L.  B.  Felsen,  "Lateral  Waves  on  an  Anisotropic  Plasma  Interface",  IRE  Trans.  PGAP 
(Communications),  AP-10 ,  No.  3,  (1962) 

L.  B.  Felsen,  "Radiation  from  a  Uniaxially  Anisotropic  Plasma  Half-Space",  to  be 
published  in  IEEE  Trans.  PGAP 

A.  Hessel,  N.  Marcuvitz  and  J.  Shmoys,  "Scattering  and  Guided  Waves  at  an  Interface 
Between  Air  and  Compressible  Plasma",  IRE  Trans.  PGAP,  AP-10,  48  - 
54,  Jan.  1962 

C.  J.  Marcinkowski,  "Some  Diffraction  Patterns  of  an  Absorbing  Half-Plane",  Appl. 
Sci.  Res.,  ^  189-198  (1961) 

C.  J.  Marcinkowski  and  L.  B.  Felsen,  "On  High-Frequency  Diffraction  by  Curved 

Surfaces  with  Periodic  Impedance  Properties",  J.  Res.  NBS,  66D,  No.  6, 
699  (1962) 

L.  B.  Felsen  and  C.  J.  Marcinkowski,  "Diffraction  by  a  Cylinder  with  a  Variable 
Surface  Impedance",  Proc.  Roy.  Soc.  London,  A267  ,  329-350  (1962) 

C.  J.  Marcinkowski  and  L.  B.  Felsen,  "On  the  Geometrical  Optics  of  Curved  Surfaces 
with  Periodic  Impedance  Properties",  J,  Res.  NBS,  66D,  699-705  (Nov. - 
Dec.  1962) 

C.  J.  Marcinkowski  and  L.  B.  Felsen,  "On  the  Limitations  of  Geometrical  Optics 
Solutions  for  Curved  Surfaces  with  Variable  Impedance  Properties", 

J.  Res.  NBS,  66D,  No.  6,  707-712  (Nov.  -  Dec.  1962) 
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Paper!  published  in  $clentiflc  jottrnal«  (continued) 

A.  A.  Oliner  and  T.  Tamir,  "Backward  Waves  on  Isotropic  Plasma  Slabs",  J.  Appl. 
Phys.,  33,  231-233,  Jan.  1962 

A.  A.  Oliner  and  T.  Tamir,  "Radiation  from  Semi-Infinite  Slot-Excited  Plasma 

Sheath  Configurations",  Second  Symposium  on  the  Plasma  Sheath,  Plenum 
Press,  to  be  published 

J.  Shmoys,  "Diffraction  by  a  Half-Plane  with  a  Special  Impedance  Variation",  IRE 
Trans.  PGAP,  AP-7,  S88-S90,  Dec.  1959 

J.  Shmoys  and  E.  Mishkin,  "Hydromagnetic  Waveguides  with  Finite  Conductivity  and 
Arbitrary  Croes-Sertion",  Phys.  and  Fluids,  2>  (I960) 

J.  Shmoys,  "A  Proposed  Diagnostic  Method  for  Cylindrical  Plasmas",  J.  Appl.  Phys., 
32,  689  (1961) 

T;  Tamir  and  A.  A.  Oliner,  "The  Influence  of  Complex  Waves  on  the  Radiation  Field 
of  a  Slot-Excited  Plasma  Layer",  IRE  Trans.  PGAP,  AP-10,  55-65,  Jan. 
1962  - 

T.  Tamir  and  A.  A.  Oliner,  "The  Spectrum  of  Electromagnetic  Waves  Guided  by  a 
Plasma  Layer",  Proc.  IEEE,  Vol.  51  317-332,  Feb.  1963 

D.  Youla,  "A  Normal  Form  for  a  Matrix  under  the  Unitary  Congruence  Group",  Can. 
Jour.  Math.,  Vol.  13,  694-704,  1961 

D.  Youla,  "On  the  Factorization  of  Rational  Matrices",  IRE  Trans,  on  Information 
Theory,  July  1961 

D.  Youla,  "New  Theory  of  Cascade  Synthesis",  IRE  Trans,  on  Circuit  Theory,  Sept. 
1961 


Talks  presented  at  technical  meetings; 

E.  Arbel  and  L.  B.  Felsen,  "Radiation  from  Sources  in  Anisotropic  Plasmas", 

Symposium  on  Electromagnetic  Theory  and  Antennas,  Copenhagen,  June 
1962 

L.  B.  Felsen,  "Electromagnetic  Properties  of  Wedge  and  Cone  Surfaces  with  a  Linearly 
Varying  Surface  Impedance",  Symposium  on  Electromagnetic  Theory,  Univ. 
of  Toronto,  July  1959 

L.  B.  Felsen,  "Diffraction  by  Objects  in  a  Certain  Variable  Medium",  URSI  Spring 
Meeting,  Washington,  D.  C. ,  1960 

L.  B.  Felsen  and  C.J.  Marcinkowski,  "Diffraction  by  a  Cylinder  with  a  Variable 
Surface  Impedance",  URSI  Spring  Meeting,  Washington,  D.  C. ,  I960 

L.  B.  Felsen,  "Lateral  Waves  on  an  Anisotropic  Plasma  Interface",  URSI  Spring 
Meeting,  Washington,  D.  C. ,  1962 
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Talkg  presented  at  technical  meetings  (continued) 

L.  B.  Felaen,  *  Two- Dimenaional  Radiation  and  Diffraction  Problems  in  Uniaxially 
Anisotropic  Media*.  URSI  Spring  Meeting.  Washington,  D.  C. ,  1962 

C.  J.  Marcinkowski,  "Plane  Wave  Scattering  by  a  Linearly  Varying  Surface  Impedance*, 
Bull.  Amer.  Phys.  Soc. ,  Series  II,  Vol.  7,  31,  Jan.  1962 

A.  A.  Oliner  and  T.  Tamir,  "Radiation  from  Semi-Infinite  Slot-Excited  Plasma  Sheath 
Configurations",  Second  Symposium  on  the  Plasma  Sheath,  Boston,  Mass. 
April  1962 

J.  Shmoys,  "Diffraction  by  a  Half-Plane  with  a  Special  Impedance  Variation", 
Symposium  on  Electromagnetic  Theory,  Toronto,  June  1959 

T.  Tamir  and  A.  A.  Oliner,  "The  Near  Field  Supported  by  a  Source-Excited  Plasma 
Slab",  URSI-IRE,  Washington,  D.  C. ,  May  1961 

T.  Tamir  and  A.  A.  Oliner,  "The  Influence  of  Complex  Waves  on  the  Radiation  Field 
of  a  Slot- Excited  Plasma  Layer",  URSI-IRE,  Austin,  Texas,  Nov.  1961 


\ 


RADC  (RAYLD)  Attn:  Document*  Library 
Griffi**  AFB.  N.  Y. 

RADC  (RCE) 

Attn:  Dr.  John  S.  Burge ■* 

Griffi**  AFB,  N,  Y. 

AF  Mi**il«  Dev.  Cent.  (MDGRT) 
Holloman  AFB,  New  Mexico 

Director  of  Resident  Training 
3380th  Technical  Training  Group 
Keesler  AFB,  Mississippi 
Attn:  OA-3011  Course 

AUL 

Maxwell  AFB,  Ala, 

AF  Missile  Test  Center 
Patrick  AFB,  Fla. 

Attn;  AFMTC,  Tech.  Library,  MU-135 


USAF  Security  Service  (CLR) 

San  Antonio,  Texas 

Commanding  General 
USASRDL 

Ft.  Monmouth,  New  Jersey 
Attn:  Tech.  Doc.  Ctr. 

SIGRA/SL-ADT 

Department  of  the  Army 
Office  of  the  Chief  Signal  Officer 
Washington  25,  D.  C. 

Attn:  SIGRD-4a-2 

Massachusetts  Institute  of  Technology 
Signal  Corps  Liaison  Officer 
Cambridge  39,  Mass. 

Attn:  A.  D.  Bedrosian,  Room  26-131 

Office  of  Chief  Signal  Officer 
Engineering  and  Technical  Division 
Washington  25,  D.  C. 

Attn;  SIGNET-5 

Director 

U.S.  Army  Ordnance 
Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland 
Attn:  Ballistic  Measurements  Laboratory 


Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland 
Attn:  Technical  Information  Branch 

Guided  Missile  Fuze  Library 
Diamond  Ordnance  Fuze  Laboratories 
Washington  25,  D.  C. 

Attn:  R.  D.  Hatcher,  Chief  Microwave 
Development  Section 

Commanding  General 
USASRDL 

Fort  Monmouth,  New  Jersey 
Attn:  SIGFM/EL-AT 

Redstone  Scientific  Information  Center 
U.S.  Army  Missile  Corhm.and 
Redstone  Arsenal,  Alabam.a  (5  copies) 

Technical  Information  Office 
European  Office,  Aerospace  Research 
Shell  Bldg. ,  47  Cantersteen 
Brussels,  Belgium 

OAR  (RROS,  Col.  John  R.  Fowler) 

Tempo  D 

4th  and  Independence  Ave. 

Washington  25,  D.  C. 

AFOSR,  OAR  (SRYP) 

Tempo  D 

4th  and  Independence  Ave. 

Washington  25,  D.  C. 

Hq.  USAF  (AFOAC-S/E) 

Comm.unications  -  Electronics  Directorate 
Washington  25,  D.  C. 

Hq.  OAR  (RROSP,  Maj.  Richard  W.  Nelson) 
Tempo  D,  4th  and  Independence  Ave. 
Washington  25,  D.  C. 

ASD  (ASAPRD  -  Diet) 

Wright-Patterson  AFB,  Ohio 

WADD  (WCLRSA,  Mr.  Fortune) 
Wright-Patterson  AFB,  Ohio 

ASD  (ASRNRE  -  3) 

Attn:  Mr.  Paul  Springer 
Wright-Patterson  Af3,  Ohio 

Director,  Electronics  Div. 

Air  Technical  Intelligence  Center 
Attn:  AFCIN-4E1,  Colonel  H.  K.  Gilbert 
Wright-Patterson  AFB,  Ohio 
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WADD  (WWDRTR,  Mr.  A.  D.  Clark) 
Directorate  of  Syetem  Engineering 
Oyna  Soar  Engineering  Office 
Wright-Pattereon  AFB,  Ohio 


National  Aeronautical  Space  Agency 
Langley  Aeronautical  Research  Laboratory 
Langley,  Virginia 
Attn;  Mr.  Cliff  Nelson 


Lt.  Col.  Jensen  (SSRTW) 

Space  Systems  Division 
Air  Force  Unit  Post  Office 
Los  Angeles  45,  California 

Director 

Evans  Signal  Laboratory 
Belmar,  New  Jersey 
Attn:  Mr.  O.  C.  Woodyard 

Commanding  General,  SIGFM/EL-PC 
USASRDL 

Fort  Monmouth,  New  Jersey 
Attn;  Dr.  Horst  H.  Kedesdy 
Deputy  Chief,  Chem-Physics  Branch 

ASTIA  (TIPA-M 

Arlington  Hall  Station  (10  copies) 

Arlington  12,  Virginia 

Library 

National  Bureau  of  Standards 

Boulder  Laboratories 

Boulder,  Colorado  (2  copies) 

Defense  Research  Member 
Canadian  Joint  Staff 
2450  Massachusetts  Avenue,  N.  W. 
Washington  8,  D.  C. 


AFCRL,  OAR  (CRXRA  -  Stop  39) 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

(10  copies) 

AFCRL,  Office  of  Aerospace  Research  (CRRD) 

Attn:  Contract  Files 

L.  G.  Hanscom.  Field,  Bedford,  Mass. 

(2  copies) 

AFCRL,  Office  of  Aerospace  Research  (CRRD) 
Attn;  Carlyle  J.  Slctten  ' 

L.  G,  Hanscon;  Field,  Bedford,  Mass. 

(3  copies) 

Hq.  ESD  (ESRDW,  Major  John  J.  Hobson) 

L.G.  Hanscom.  Field 
Bedford,  Mass. 

Electronic  Systems  Division  (AFSC) 

Technical  Information  Services  Div.  (ESAT) 
L.G.  Hanscom,  Field,  Bedford,  Mass. 

Hq.  AFCRL.  OAR  (CRXR,  J.  R.  Marple) 

L.G.  Hanscom  Field,  Bedford,  Mass. 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 
Attn:  Code  690 


National  Bureau  of  Standards 
U.  S.  Departm.ent  of  Comm.erce 
Washington  25,  D.  C. 

Attn;  Mr,  A.  G,  McNish 


Chief,  Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  DLI-31  (2  copies) 


National  Bureau  of  Standards 
U.  S.  Department  of  Commerce 
Washington  25,  D.  C. 

Attn;  Gustave  Shapiro,  Chief 
Engineering  Electronics  Section 
Electricity  and  Electronics  Div. 

Office  of  Scientific  Intelligence 
Central  Intelligence  Agency 
2430  E  Street,  N.  W. 

Washington  25,  D.  C. 

Director 

National  Security  Agency 

Fort  George  G.  Meade,  Maryland 

Attn:  C3/TDL 


Com.mander 

U.S.  Naval  Air  Missile  Test  Center 
Point  Mugu,  California 
Attn;  Code  366 

U.S.  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring  19,  Maryland 
Attn:  The  Library 

Commander 

U.S.  Naval  Ordnance  Test  Station 
China  Lake,  California 
Attn:  Code  753 
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librarian 

U.  S.  Naval  Postgraduate  School 
Monterey,  California 

National  Aeronautics  and  Space  Administration 


Attn:  Antenna  Systems  Branch 
Goddard  Space  Flight  Center 
Greenbelt,  Maryland 

Director 

U.S.  Naval  Research  Laboratory 
Washington  25,  O.  C. 

Attn;  Code  2027  (  2  copies  ) 

Dr.  J.  I.  Bohnert,  Code  5210 
U.S.  Naval  Research  Laboratory 
Washington  25,  D.  C. 


AFSC  Scientific  and  Technical  Liaison 
Office 

c/o  Department  of  the  Navy 
Room  2  305,  Munitions  Bldg. 
Washington  25,  D.  C. 

Aero  Geo  Astro  Corp. 

1200  Duke  Street 
Alexandria,  Virginia 
Attn:  Library 


Commanding  Officer  and  Director 

U.  S.  Navy  Underwater  Sound  Laboratory 

Fort  Trumbull,  New  London,  Conn. 

Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  427 

Commanding  Officer 

U.S.  Naval  Air  Development  Center 

Johnsville,  Penn. 

Attn;  NADC  Library 

Office  of  Naval  Research 
Branch  Office,  London 
Navy  100,  Box  39 

F.P.O.  New  York,  N.  Y.  (10  copies) 

Commanding  Officer  and  Director 

U.S.  Navy  Electronics  Laboratory  (Library) 

San  Diego  52,  Calif. 

Commander 

U.S.  Naval  Air  Test  Center 
Patuxent  River,  Maryland 
Attn:  ET-315,  Antenna  Branch 

Commanding  Officer 

U.S.  Naval  Ordnance  Laboratory 

Corona,  Calif. 

Attn:  Documents  Librarian 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  817B 


Aerospace  Corp. 

Box  95085 

Los  Angeles  45,  Calif. 

Attn:  Library 

Airborne  Instruments  Laboratory,  Inc. 
Division  of  Cutler  Hammer 
Walt  Whitman  Road 
Melville,  L.  I. ,  N.  Y. 

Attn;  Library 

Aircom,  Inc. 

48  Cum,mington  Street 
boston.  Mass. 

Andrew  Alford,  Consulting  Engineers 
299  Atlantic  Ave. 

Boston  10,  Mass. 

Aerospace  Corp. 

Satellite  Control 
Attn;  Mr.  R.  C.  Hansen 
Post  Office  Box  95085 
Los  Angeles  45,  Calif. 

ACF  Electronics  Division 
Bladensburg  Plant 
52nd  Avenue  and  Jackson  Street 
Bladensburg,  Maryland 
Attn;  Librarian 

Battelle  Memorial  Institute 
505  King  Ave. 

Columbus  1,  Ohio 

Attn:  Wayne  E.  Rife,  Project  Leader 
Electrical  Engineering  Division 
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Bell  Aircraft  Corporation 
Poet  Office  Box  One 
Buffalo  5,  N.  Y. 

Attn:  Eunice  P.  Hazelton,  Librarian 

Bell  Telephone  Laboratoriee 
Murray  Hill 
New  Jersey 

Bell  Telephone  Laboratories,  Inc. 
Technical  Information  Library 
Whippany  Laboratory 
Whippany,  New  Jersey 
Attn:  Technical  Reports  Librarian 

Beiidix  Pacific  Division 
11600  Sherman  Way 
North  Hollywood,  Calif. 

Attn:  Engineering  Library 

Bendlx  Radio  Division 
Bendix  Aviation  Corporation 
E.  Joppa  Road 
Towson  4,  Maryland 
Attn:  Dr.  D,  M,  Allison,  Jr. 

Director  Engineering  and  Research 

BJorksten  Research  Laboratories,  Inc. 

P.  O.  Box  265 
Madison,  Wisconsin 
Attn:  Librarian 

Boeing  Airplane  Company 

Pilotless  Aircraft  Division 

P.  O.  Box  3707  (2  copies) 

Seattle  24,  Washington 

Attn:  R.  R.  Barber,  Library  Supervisor 

Boeing  Company 
3801  S.  Oliver  Street 
Wichita  1,  Kansas 

Attn:  Kenneth  C,  Knight,  Library  Supervis 

Chance  Vought  Corp. 

9314  West  Jefferson  Blvd. 

Dallas,  Texas 

Attn;  A.  D.  Pattullo,  Librarian 

Chance  Vought  Corporation 
Vought  Electronics  Div. 

P.  O.  Box  5907 
Dallas  22,  Texas 


Chu  Associates 
P.  O.  Box  387 
Whitcomb  Avenue 
Littleton,  Mass. 

Collins  Radio  Co. 

855  35th  Street,  N.  E. 

Cedar  Rapids,  Iowa 
Attn;  Dr.  R.  L.  McCreary 

Convair,  A  Div.  of  General  Dynamics  Corp. 
Fort  Worth,  Texas 
Attn;  K.  G.  Brown 

Division  Research  Librarian 

Convair,  A,  Div,  of  General  Dynamics  Corp. 
3165  Pacific  Hwy. 

San  Diego  12,  California 
Attn;  Mrs.  Dora  B.  Burke 
Engineering  Librarian 

Cornell  Aeronautical  Laboratory  Inc. 

4455  Genesee  Street 
Buffalo  21,  New  York 
Attn;  Librarian 

Dalmo  Victor  Company 
A  Division  of  Textron,  Inc. 

1515  Industrial  Way 
Belmont,  Calif. 

Attn;  Mary  Ellen  Addems,  Technical 
Librarian 

Dome  and  Margolin,  Inc. 

29  New  York  Avenue 
Westbury,  L.  I. ,  N.  Y. 


or 

Aircraft  Division 
Douglas  Aircraft  Co. ,  Inc. 

3855  Lakewood  Blvd. 

Long  Beach,  California 
Attn;  Technical  Library 

Douglas  Aircraft  Co. ,  Inc. 

3000  Ocean  Park  Blvd. 

Santa  Monica,  California 
Attn:  Peter  Duyan,  Jr. 

Chief,  Electrical/Electronics  Section 
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Douglai  Aircraft  Co. ,  Inc. 

2000  North  Memorial  Drive 
Tulsa,  Oklahoma 

Attn;  Engineering  librarian,  0-250 

Electromagnetic  Research  Corporation 
5001  College  Ave. 

College  Park,  Maryland 
Attn:  Mr.  Martin  Katzin 

Electronics  Communication 
1830  York  Road 
Timonium,  Maryland 

Electronic  Specialty  Company 
5121  San  Fernando  Road 
Los  Angeles  39,  Calif. 

Attn:  Donald  L.  Margerum 

Chief  Engineer,  Radiating  Systems 
Div, 

Emerson  and  Cuming,  Inc. 

59  Walpole  Street 
Canton,  Maas. 

Attn:  Mr.  W.  Cuming. 

Emerson  Electric  Mfg.  Co. 

8100  West  Florissant  Avenue 

St.  Louis  21,  Missouri 

Attn:  Mr.  E.  R.  Breslin,  Librarian 

Fairchild  Aircraft  -  Missiles  Division 
Fairchild  Eng.  and  Airplane  Corp. 
Hagerstown,  Maryland 
Attn:  Library 

ITT  Federal  Laboratories 
Technical  Library 
500  Washington  Avenue 
Nutley  10,  New  Jersey 

Gabriel  Electronics  Division 
Main  and  Pleasant  Streets 
Millie,  Mass. 

Attn:  Dr.  Edward  Altshuler 


General  Electric  Company 
Missile  and  Space  Vehicle  Dept. 

3198  Chestnut  St. ,  Philadelphia,  Penn. 
Attn;  Documents  Library 

General  Electric  Company 
3750  D  Street 
Philadelphia  24,  Pa. 

Attn:  Mr.  H.  G.  Lew 

Missile  and  Space  Vehicle  Dept. 

General  Precision  Laboratory,  Inc. 

63  Bedford  Road 
Pleasantville,  New  York 
Attn:  Librarian 

Goodyear  Aircraft  Corp. 

1210  Massillon  Road 
Akron  15,  Ohio 
Attn:  Library,  Plant  G 

Granger  Associates 
Electronic  Systems 
974  Commercial  Street 
Palo  Alto,  Calif. 

Attn:  John  V.  N.  Granger,  President 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,  L.  I. ,  N.  Y. 

Attn;  Engineering  Librarian,  Plant  No.  5 

Hallic  rafters  Company 
4401  West  5th  Ave, 

Chicago  24,  Illinois 

Attn:  La  Verne  LaGioia,  Librarian 

The  Hallicrafters  Co, 

5th  and  Kostner  Avenues 
Chicago  24,  Illinois 
Attn:  Henri  Hodara,  Head  of  Space 
Communication 

Hoffman  Electronics  Corp. 

3761  South  Hill  Street 
Los  Angeles  7,  Calif. 

Attn:  Engineering  Library 


General  Electric  Company 
Electronics  Park 
Syracuse,  New  York 
Attn;  Documents  Library 

B.  Fletcher,  Bldg.  3-143A 


Hughes  Aircraft  Company 
Antenna  Department 
Bldg.  12,  Mail  Station  2714 
Culver  City,  Calif. 

Attn:  Dr.  W.  H.  Kummer 
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Hughea  Aircraft  Company 
Florence  Ave.  and  Teale  Street 
Culver  City.  Calif. 

Attn:  Louis  L.  Bailin 

Manager,  Antenna  Dept. 

Hughea  Aircraft  Company 

Attn:  Mr.  L.  Stark,  Microwave  Dept. 

Radar  Laboratory,  P.  O  Box  2097 
Bldg.  600,  Mail  Station  C-152 
Fullerton,  Calif. 

International  Business  Machines  Corp. 

Space  Guidance  Center  -  Federal  Systems  Div. 
Owego,  Tioga  County,  New  York 
Attn:  Technical  Reports  Center 

International  Resistance  Company 
401  N.  Broad  Street 
Philadelphia  8,  Pa. 

Attn:  Research  Library 

ITT  Federal  Laboratories 
3700  East  Pontiac  Street 
Fort  Wayne  1,  Indiana 
Attn:  Technical  Library 

Atlantic  Research  Corporation 
Shirley  Highway  at  Edsall  Road 
Alexandria,  Virginia 
Attn:  Delmer  C.  Ports 

Dr.  Henry  Jasik,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Lockheed  Aircraft  Corporation 
2555  N.  Hollywood  Way 
California  Division  Engineering  Library 
Department  72-25,  Plant  A-1,  Bldg.  63-1 
Burbank,  Calif. 

Attn:  N.  C.  Harnois 

Lockheed  Aircraft  Corporation 
Missiles  and  Space  Division 
Technical  Information  Center 
3251  Hanover  Street 

Palo  Alto,  Calif.  (2  copies) 

Martin-Marietta  Corp. 

12250  S.  State  Highway  65, 

Jefferson  County,  Colorado 
Attn:  Mr.  Jack  McCormick 


The  Martin  Co. 

Baltimore  3,  Maryland 
Attn;  Engineering  Library 
Antenna  Design  Group 

Mathematical  Reviews 
l90  Hope  Street 
Providence  6,  Rhode  Island 

The  W.  L.  Maxson  Corporation 
475  10th  Avenue 
New  York,  N.  Y. 

Attn:  Miss  Dorothy  Clark 

McDonnell  Aircraft  Corp. ,  Dept.  644 
Box  516,  St.  Louis  66,  Missouri 
Attn:  C.  E.  Zoller  ' 

Engineering  Library 

McMillan  Laboratory,  Inc. 

Brownyille  Avenue 
Ipswich,  Mass. 

Attn;  Security  Officer,  Document  Room 

Melpar,  Inc. 

3000  Arlington  Blvd. 

Falls  Church,  Virginia 

Attn;  Engineering  Technical  Library 

Microwave  Associates,  Inc. 

South  Avenue 
Burlington,  Mass. 

Microwave  Development  Laboratories, 
Inc. 

92  Broad  Street 
Wellesley  57,  Mass. 

Attn;  N.  Tucker,  General  Manager 

The  Mitre  Corporation 
244  Wood  Street 
Lexington  73,  Mass. 

Attn;  Mrs.  Jean  E.  Claflin,  Librarian 

Motorola,  Inc. 

8201  East  McDowell  Road 

Phoenix,  Arizona 

Attn:  Dr.  Thomas  E.  Tice 

Motorola,  Inc. 

Phoenix  Research  Laboratory 
3102  N.  56th  Street 
Phoenix,  Arizona 
Attn:  Dr.  A.  L.  Aden 
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National  Research  Council 
Radio  and  Electrical  Engineering  Div. 
Ottawa,  Ontario,  Canada 
Attn:  Dr,  G.A.  Miller,  Head 
Microwave  Section 

North  American  Aviation,  Inc. 

12214  Lakewood  Blvd. 

Downey,  Calif. 

Attn:  Technical  Information  Center  (495-12) 
Space  and  Information  Systems  Div. 

North  American  Aviation,  Inc. 

Los  Angeles  International  Airport 
Los  Angeles  45,  Calif. 

Attn:  Engineering  Technical  File 

Page  Communications  Engineers,  Inc. 

2001  Wisconsin  Avenue,  N,  W. 

Washington  7,  D,  C. 

Attn:  (Mrs.)  Ruth  Temple,  Librarian 

Northrop  Corporation 
Norair  Division 
1001  East  Broadway 
Hawthorne,  Calif, 

Attn:  Technical  Information  3924-31 

Philco  Corporation 
Research  Division 
Union  Meeting  Pond 
Blue  Bell,  Pa, 

Attn:  Research  Librarian 

Pickard  and  Burns,  Inc, 

103  Fourth  Avenue 
Waltham  54,  Mass. 

Attn:  Dr.  Richard  H.  Woodward 

Polytechnic  Research  and  Development  Co. 
202  Tillary  Street 
Brooklyn  1,  N.  Y. 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass, 

Attn:  Dr.  John  Ruze 

Radiation,  Inc. 

Melbourne,  Florida 
Attn:  RF  Systems  Div. 

Technical  Information  Center 


Radiation  Systems,  Inc. 

440  Swann  Avenue 
Alexandria,  Virginia 
Attn:  Library 

RCA  Laboratories 

David  Sarnoff  Research  Center 

201  Washington  Road 

Princeton,  New  Jersey 

Attn:  Miss  Fern  Cloak,  Librarian 

Radio  Corp.  of  America 
Defense  Electronic  Products 
Bldg.  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader,  Staff 
Engineer 

Organization  of  Chief  Technical 
Adm.inistrator 

Radio  Corp,  of  America 
Missile  Control  and  Electronics  Div. 
Bedford  Street,  Burlington,  Mass. 
Attn:  Librarian 

Radio  Corp.  of  America 
West  Coast  Missile  and  Surface  Radar 
Div. 

Engineering  Library,  Bldg.  306/2 

Attn:  L.  R.  Hund,  Librarian 

8500  Balboa  Blvd.  ,  Van  Nuys,  Calif. 

Radio  Corp,  of  America 
Defense  Electronic  Products 
Advanced  Military  Systems 
Princeton,  Now  Jersey 
Attn:  Mr.  David  Shore 

Director,  USAF  Project  RAND 
Via:  AF  Liaison  Office 
The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  Calif. 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  Calif. 

Attn:  Technical  Library 

Rantec  Corporation 
23999  Ventura  Blvd. 

Calabasas,  Calif. 

Attn:  Grace  Keener,  Office  Manager 
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Raytheon  Company 
Boaton  Poet  Road 
Wayland,  Mate. 

Attn;  Mr.  Robert  Borte 

Raytheon  Company 
Wayland  Laboratory 
Wayland,  Maas. 

Attn:  Miss  Alice  G.  Anderson,  Librarian 

Raytheon  Company 
Missile  Systems  Div. 

Hartwell  Road,  Bedford,  Mass. 

Attn:  Donald  H.  Archer 

Remington  Rand  UNIVAC 
Div.  of  Sperry  Rand  Corp. 

P.  O.  Box  500 
Blue  Bell,  Pa. 

Attn:  Engineering  Library 

Republic  Aviation  Corp. 

Farmingdale,  L.  I. ,  N.  Y. 

Attn:  Engineering  Library 

Thru:  AF  Plant  Repr.  Republic 
Aviation  Corp. 
Farmingdale,  L.  I. ,  N.  Y. 

Ryan  Aeronautical  Company 
2701  Harbor  Drive 
Lindbergh  Field 
San  Diego  12,  Calif. 

Attn:  Library 

Sage  Laboratories,  Inc. 

3  Huron  Drive 
Natick,  Mass. 

Sanders  Associates,  Inc. 

95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  Mr.  Norman  R.  Wild 

Sandia  Corporation 
P.O.  Box  5800 
Albuquerque,  New  Mexico 
Attn:  Records  Management  and 
Services  Department 

Scanwell  Laboratories,  Inc. 

6601  Scanwell  Lane 
Springfield,  Va. 


STL  Technical  Library 

Document  Acquisitions 

Space  Technology  Laboratories,  Inc. 

P.O.  Box  95001 

Los  Angeles  45,  California 

Sperry  Gyroscope  Company 
Great  Neck,  L.  I. ,  N.  Y. 

Attn;  Florence  W.  Turnbull 
Engineering  Librarian 

Stanford  Research  Institute 
Documents  Center 
Menlo  Park,  Calif. 

Attn;  Acquisitions 

Sylvania  Electric  Products,  Inc. 

100  First  Ave.,  Waltham  54,  Mass. 
Attn:  Charles  A.  Thornhill 
Report  Librarian 
Waltham  Laboratories  Library 

Sylvania  Elec,  Prod.  Inc. 

Electronic  Defense  Laboratory 
P.O.  Box  205 
Mountain  View,  Calif. 

Attn:  Library 

Sylvania  Reconnaissance  Systems  Lab. 
Box  188,  Mountain  View,  Calif. 

Attn:  Marvin  D.  Waldman 

TRG,  Inc. 

400  Border  Street,  East  Boston,  Mass. 
Attn:  Dr.  Alan  F.  Kay 

A.  S,  Thomas,  Inc. 

355  Providence  Highway 
We3t.vood,  Mass. 

Attn:  A.  S.  Thomas,  President 

Texas  Instruments,  Inc. 

6000  Lemmon  Avenue 
Dallas  9,  Texas 
Attn;  John  B.  Travis 

Systems  Planning  Branch 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l.  Airport  Box  1897 
Baltimore  3,  Maryland 
Attn;  Engineering  Library 
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Library  Geophysical  Institute  of  the 
University  <A  Alaska 
College,  Alaska 

Brown  University 

Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Or.  C.  M.  Angulo 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  Calif. 

Attn:  Mr.  I.  E.  Newlan 

California  Institute  of  Technology 
1201  E.  California  Drive 
Pasadena,  Calif. 

Attn;  Dr.  C.  Papas 

Space  Sciences  Laboratory 
Leuschner  Observatory 
University  of  California 
Berkeley  4,  Calif. 

Attn;  Or.  Samuel  Silver,  Prof,  of  Engrg. 
Science  and  Director,  Space 
Sciences  Laboratory 

University  of  California 
Electronics  Research  Lab. 

332  Cory  Hall,  Berkeley  4,  Calif. 

Attn;  J.  R.  Whinnery 

University  of  Southern  California 
University  Park 
Los  Angeles,  California 
Attn;  Dr.  Raymond  L.  Chuan 

Director,  Engineering  Center 

Case  Institute  of  Technology 
Electrical  Engineering  Dept. 

10900  Euclid  Avenue 

Cleveland,  Ohio 

Attn;  Professor  Robert  Plonsey 

Columbia  University 
Dept,  of  Electrical  Engineering 
Morningside  Heights,  New  York,  N.  Y. 
Attn;  Dr.  Schlesinger 

University  of  Southern  California 
University  Park 
Los  Angeles  7,  Calif. 

Attn;  Z.A.  Kaprielian,  Assoc.  Prof,  of 
Electrical  Engineering 


Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Dalman 

University  of  Florida 

Dept,  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:  Prof.  M.  H.  Latour,  Library 

Library 

Georgia  Technology  Research  Institute 
Engineering  Experiment  Station 
722  Cherry  Street,  N.  V, 

Atlanta,  Georgia 

Attn:  Mrs.  J.  H.  Crosland,  Xibrarian 

Harvard  University 
Technical  Reports  Collection 
Gordon  McKay  Library 
303  Pierce  Hall 

Oxford  Street,  Cambridge  38,  Mass. 

Attn:  Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,  Mass. 

Attn:  Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,  Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,  Illinois 

Attn:  Dr.  P,  E.  Mayes,  Dept,  of  Electrical 
Engineering 

Illinois  Institute  of  Technology 

Technology  Center 

Dept,  of  Electrical  Engineering 

Chicago  16,  Illinois 

Attn:  Paul  C.  Yuen 

Electronics  Research  Laboratory 

The  John  Hopkins  University 

Homewood  Campus 

Baltimore  18,  Maryland 

Attn;  Dr.  Donald  E.  Kerr,  Dept,  (vf  Physics 

The  John  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue,  Silver  Spring,  Md. 
Attn;  Mr.  George  L.  Seielstad 
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University  of  Kansas 
Electrical  Engineering  Department 
Lawrence,  Kansas 
Attn:  Dr.  H.  Unz 

Ijowell  Technological  Institute 

Research  Foundation 

P.  O.  Box  709i  Lowell,  Mass. 

Attn:  Dr.  Charles  R.  Mingins 

Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Bldg.  26,  Room  327 
Cambridge  39,  Mass. 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology 
Lincoln  Laboratory 
P.O.  Box  7  3 
Lexington  73,  Mass. 

Attn:  Mary  A.  Granese,  Librarian 

McGill  University 

Dept,  of  Electrical  Engineering 

Montreal,  Canada 

Attn:  Dr.  T.  Pavlasek 

University  of  Michigan 
Electronic  Defense  Group 
Institute  of  Science  and  Technology 
Ann  Arbor,  Michigan 
Attn:  J.  A.  Boyd,  Supervisor 

University  of  Michigan 

Office  of  Research  Administration 

Radiation  Laboratory 

912  N.  Main  Street 

Ann  Arbor,  Michigan 

Attn:  Mr.  Ralph  E.  Hiatt 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories,  Willow  Run 
Ypsilanti,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  14,  Minnesota 

Attn:  Mr.  Robert  H.  Stum.m,  Library 

Physical  Science  Laboratory 
New  Mexico  State  University 
University  Park,  New  Mexico 
Attn:  Mr.  H.  W.  Haas 


Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E.  Beam 

Antenna  Laboratory 

Dept,  of  Electrical  Engineering 

The  Ohio  State  University 

2024  Neil  Avenue 

Columbus  10,  Ohio 

Attn.  Reports  Librarian 

The  University  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 
Attn:  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

University  of  Pennsylvania 
Institute  of  Cooperative  Research 
3400  Walnut  Street,  Philadelphia,  Pa. 
Attn;  Dept,  of  Electrical  Engineering 

Polytechnic  Ins'  '  >fe  of  Brooklyn 
Microwave  Resi\i  h  Institute 
55  Johnson  Street,  Brooklyn,  N.  Y. 
Attn:  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn,  N.  Y. 

Attn:  Mr.  A.  E.  Laem.mel 

The  Pennsylvania  State  University 
223  Electrical  Engineering 
University  Park,  Pa. 

Attn:  A.  H.  Waynick,  Director 
Ionosphere  Research  Lab. 


Airport 

Purdue  University 

Dept,  of  Electrical  Engineering 

Lafayette,  Ind. 

Attn;  Dr.  Schultz 

Library 

W.  W.  Hansen  Laboratory  of  Physics 
Stanford  University 
Stanford,  Calif. 
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Syracuae  Unlvaraity  Reaeareh  Institute 
Collendale  Campus 
Syracuse  10,  N.  Y. 

Attn:  Dr.  C.S.  Grove.  Jr. 

Director  of  Engineering  Research 

Technical  University 
destervoldgade  10  G 
Cq;>enhagen,  Denmark 
Attn:  Prof.  Hans  Lottrup  Knudsen 

University  of  Tennessee 
Ferris  Hall 

W.  Cumberland  Avenue 
Knoxville  16,  Tennessee 

The  University  of  Texas 
Electrical  Engineering  Research  Lab. 

P.  O.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhardt 
Assistant  Director 

The  University  of  Texas 
Defense  Research  Laboratory 
Austin,  Texas 
Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 
Dept,  of  Electrical  Engineering 
Toronto,  Canada 
Attn;  Prof.  G.  Sinclair 

University  of  Washington 
Dept,  of  EJlectrical  Engineering 
Seattle  5,  Washington 
Attn;  D.  K.  Reynolds 

University  of  Wisconsin 
Dept,  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn;  Dr.  Scheibe 

WADD  (WCLKTR) 

Wright-Patterson  AFB,  Ohio 

U.  S.  Atomic  Energy  Commission 
Office  of  Technical  Information  Extension 
P.  O.  Box  62 
Oak  Ridge,  Tennessee 


Harvard  University 
Cordon  McKay  Laboratory 
9  Oxford  Street 
Cambridge  38,  Maes, 

Attn:  Ronald  W.  P.  King 

Professor  of  Applied  Physics 

Harvard  University 
Gordon  McKay  Laboratory 
9  Oxford  Street 
Cambridge  38,  Mass. 

Attn:  Prof.  S.  R.  Seshadri 
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